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Abstract. The notion of viscosity solutions of scalar fully nonlinear partial differ- 
ential equations of second order provides a framework in which startling comparison 
and uniqueness theorems, existence theorems, and theorems about continuous de- 
pendence may now be proved by very efficient and striking arguments. The range of 
important applications of these results is enormous. This article is a self-contained 
exposition of the basic theory of viscosity solutions. 



Introduction 

The theory of viscosity solutions apphes to certain partial differential equations 
of the form F{x, u, Du, D'^u) = where : x M x x S{N) R and S{N) is 
the set of symmetric N x N matrices. The primary virtues of this theory are that 
it allows merely continuous functions to be solutions of fully nonlinear equations 
of second order, that it provides very general existence and uniqueness theorems 
and that it yields precise formulations of general boundary conditions. Moreover, 
these features go hand-in-hand with a great flexibility in passing to limits in various 
settings and relatively simple proofs. In the expression F(x, u, Du, D^u) u will be 
a real- valued function defined on some subset O of R^, Du corresponds to the 
gradient of u and D'^u corresponds to the matrix of second derivatives of u. How- 
ever, as explained below, Du and D^u will not have classical meanings and, in fact, 
the examples exhibited in §1 will convince the reader that the theory encompasses 
classes of equations that have no solutions that are differentiable in the classical 
sense. 
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In order that the theory apply to a given equation F = 0, we will require F to 
satisfy a fundamental monotonicity condition 

(0.1) F{x,r,p,X) < F{x,s,p,Y) whenever r < s and F < X; 

where r, s € K, a;,p e K^, X,Y & <5(-/V) and <S(A'') is equipped with its usual order. 
Regarding (0.1) as made up of the two conditions 

(0.2) F{x,r,p,X) < F{x,s,p,X) whenever r < s, 

and 

(0.3) F{x,r,p,X) <F{x,r,p,Y) whenever F < X, 

we will give the name "degenerate ellipticity" to the second. That is, F is said 
to be degenerate elliptic if (0.3) holds. When (0.2) also holds (equivalently, (0.1) 
holds), we will say that F is proper. 

The examples given in §1 will illustrate the fact that the antimonotonicity in 
X is indeed an "ellipticity" condition. The possibility of "degeneracies" is clearly 
exhibited by considering the case in which F{x,r,p,X) does not depend on X — 
it is then degenerate elliptic. The monotonicity in r, while easier to understand, 
is a slightly subtle selection criterion that, in particular, excludes the use of the 
viscosity theory for first order equations of the form b{u)ux = f{x) in K when b is 
not a constant function, since then F{x,r,p) = b{r)p — f{x) is not nondecreasing 
in r for all choices of p (scalar conservation laws are outside of the scope of this 
theory). 

The presentation begins with §1, which, as already mentioned, provides a list of 
examples. This rather long list is offered to meet several objectives. First, we seek 
to bring the reader to our conviction that the scope of the theory is quite broad while 
providing a spectrum of meaningful applications and, at the same time, generating 
some insight as regards the fundamental structural assumption (0.1). Finally, in 
the presentation of examples involving famous second order equations, the very 
act of writing the equations in a form compatible with the theory will induce an 
interesting modification of the classical viewpoint concerning them. 

In §2 we begin an introductory presentation of the basic facts of the theory. The 
style is initially leisurely and expository and technicalities are minimized, although 
complete discussions of various key points are given and some simple arguments in- 
conveniently scattered in the literature arc presented. Results arc illustrated with 
simple examples making clear their general nature. Section 2 presents the basic no- 
tions of solution used in the theory, the analytical heart of which lies in comparison 
results. Accordingly, §3 is devoted to explaining comparison results in the simple 
setting of the Dirichlet problem; roughly speaking, they are proved by methods 
involving extensions of the maximum principle to semicontinuous functions. Once 
these comparison results arc established, existence assertions can be established by 
Perron's method, a rather striking tale that is told in §4. With this background 
in hand, the reader will have an almost complete (sub)story and with some effort 
(but not too much!) should be able to absorb in an efficient way some of the more 
technical features of the theory that are outlined in the rest of the paper. 

Other important ideas are to be found in §6, which is concerned with the is- 
sue of taking limits of viscosity solutions and applications of this and in §7, which 
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describes the adaptation of the theory to accommodate problems with other bound- 
ary conditions and problems in which the boundary condition cannot be strictly 
satisfied. In the later case, the entire problem has a generalized interpretation for 
which there is often existence and uniqueness. While the description of these re- 
sults is deferred to §7, they are fundamental and dramatic. For example, if G{p, X) 
is uniformly continuous, degenerate elliptic, and independent of x and O C is 
open and bounded, n{x) is the unit exterior normal at the point x in its smooth 
boundary dQ and / G C{Q), then the Neumann problem 

u + G{Du,D^u) - f{x) =0 mCt, u„ = on 90 

has a unique properly interpreted solution (which may not satisfy u„ = on dfl). 

Sections 5, 8, 9 discuss variations of the basic material and need not be read in 
sequence. Section 10 is devoted to a commentary about applications (which are not 
treated in the main text), references, and possible lines of future development of the 
subject. We conclude with an appendix where the reader will find a self-contained 
presentation of the proof of the analytical heart on the presentation we have chosen. 

References are not given in the main text (with the exception of §10), but are 
to be found at the end of each section. In particular, the reader should look to the 
end of a section for further comments, references that contain details ommitted in 
the main text, and technical generality, historical comments, etc. The references, 
while numerous, are not intended to be complete, except that we have sought 
to represent all the major directions of research and areas of application. There 
are original aspects of the current presentation and the reader will note differences 
between the flavor and clarity of our presentation and that of many of the citations. 
However, equipped with the view presented here, we hope and expect that perusing 
the amazing literature that has so quickly matured will be a much more rewarding 
and efficient endeavor. 

We are grateful to R. Dorroh, M. Kocan and A. Swiech for their kind help in 
reducing the number of errors herein. 
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1. Examples 

We will record here many examples of degenerate elliptic equations mentioning, 
when appropriate, areas in which they arise. The reader is invited to scan the list 
and pause where interested it is possible to proceed to §2 at any stage. Below we 
will say either that a function F{x, r,p, X) is degenerate elliptic or that the equation 
F{x,u,Du,D'^u) = is degenerate elliptic or that the "operator" or expression 
F{x, u, Du, D'^u) is degenerate elliptic and always mean the same thing, i.e., (0.3) 
holds; the term "proper" is used in a similar fashion. 

Example 1.1. Laplace's equation. We revisit an old friend, the equation 

(1.1) -Au + c(x)ti = /(a;) 

— note the sign in front of the Laplacian. The corresponding F is naturally given 
by F(x, r,p, X) — — trace(X) + c{x)r — f(x), which is proper if c > 0. 

Example 1.2. Degenerate elliptic linear equations. Example 1.1 immediately ex- 
tends to the more general linear equation 

(1.2) - ai,j{x)-g^ + J2''i{x)-g^+c{x)uix)=f{x) 

i,j=l * ^ i=l ' 

where the matrix A(x) = {ai,j{x)} is symmetric; the corresponding F is 

N 

(1.3) F{x,r,p,X) = - tTace{A{x)X) +'Ybi{x)pi + c{x)r - f{x). 

i=l 

In this case, F is degenerate elliptic if and only if A(x) > and it is proper if also 
c{x) > 0. In the event that there are constants A, A > such that XI < A{x) < AI 
for all x where I is the identity matrix, F is said to be uniformly elliptic. 

Of course, the linear equation in divergence form 

- ^ ("^'^^''^S^) + ^^^^''^fe + cix)u{x) = fix) 



can be written as above with 

N / N 

dx. 



F{x,r,p,X) = - ti&ce{A{x)X) + ^ ( bj{x) - ^ ^(x) ] pj + c{x)r - f{x) 

provided that the indicated derivatives of the exist. 

We leave the interesting class of linear equations to turn to the totally degenerate 
case of first order equations. 



USER'S GUIDE TO VISCOSITY SOLUTIONS 



5 



Example 1.3. First order equations. The main point is that a first order oper- 
ator F{x,u,Du) is always degenerate elliptic and thus it is proper if and only if 
F(x, r,p) is nondecreasing in r G M. Proper equations of the form F{x, u, Du) = 
play a fundamental role in the classical Calculus of Variations and in Optimal Con- 
trol Theory of ordinary differential equations; in this context they are often called 
Bellman or Hamilton- Jacobi equations and then F(x, r,p) is convex in (r,p). These 
equations, in the full generality of nearly arbitrary proper functions F, are also 
crucial in Differential Games Theory where they are known as Isaacs's equations. 

Example 1.4. Quasilinear elliptic equations in divergence form. The usual notion 
of ellipticity for equations of the form 

JV Q 

(1.4) -^ — {ai{x,Du)) + b{x,u,Du)=Q 

i=i 

is the monotonicity of the vector field a{x,p) in p & mapping from to K.^. If 
enough regularity is available to carry out the differentiation, we write (1.4) as 

JV „ „2 ^ r) 

(") -E |j(«.«»)^+''(«.».««)'i:£(»=.^«)=« 

and correspondingly set 

(1.6) F{x,r,p,X) = - tTace{{Dpa{x,p))X) + b{x,r,p) - ^ ^{x,p). 

1=1 * 

The monotonicity of a in p is precisely the condition that guarantees that F is 
degenerate elliptic, and then it is proper provided we ask that b be nondecreasing 
in r. 

Two well-known instances are provided by the equation of minimal surfaces in 

nonparamctric form and the "m-Laplacc's" equation that are given, respectively, 
by a{x,p) = (1 + |pp)~^/2p, b = b{x,u) and a{x,p) = \p\"^~^p, b — b{x,r) where 
m S (1, oo). Computations show that the corresponding F's are, respectively, 

F{x,r,p,X) = -{l + \p\'^)~^ trace(X) + (1 -h |p|2)"^ trace((p(g)p)X) + 6(a;,r) 

and 

F{x,r,p,X) = -1^1™-^ trace(X) - (771-2)1^1"-* tT&ce{{p^p)X) +b{x,r). 
Example 1.5. Quasilinear elliptic equations in nondivergence form. The equation 

(1.7) - ^ a^j{x,p)-^-^ +b{x,u,Du) = 0, 

where A{x,p) = {aij{x,p)} £ S{N), contains all of the above as special cases and 
corresponds to 

F{x,r,p,X) = — tvace{A{x,p)X) + b{x,r,p), 
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which is proper ii A > and b is nondecreasing with respect to r. Two relevant 
special cases are 

—uAu + f{x, u, Du) = 

with V > and / nondecreasing in m, which may be regarded as a first-order 
Hamilton- Jacobi equation perturbed by an additional "viscosity" term —v/S.u (equa- 
tions of this type arise in optimal stochastic control), and the Levi's equation 

[dxl^ dxl) y-^\dxs) ) dxl[\dxj ^[8x2) J 

d'^u f du (hi du \ d^u f du du du 
+ 2-——[—- -—\+2 



dxidxs \dx3 dxi 8x2 J dx2dxz \dxzdx2 dxi 

which is the nonparametric formulation for a hypersurface in with vanishing 
Levi's form. Note that in this example F = — trace(A(p)X) where 

Mp) = 1+pI P3P2 + Pi 

\ P3P1 - P2 P3P2 +Pl pI+pI 

so that A > but det(A(p)) = for all p. 

Example 1.6. Hamilton- Jacobi-Bellman and Isaacs equations. Hamilton- Jacobi- 
Bellman and Isaacs equations are, respectively, the fundamental partial differential 
equations for stochastic control and stochastic differential games. The natural 
setting involves a collection of elliptic operators of second-order depending either 
on one parameter a (in the Hamilton-Jacobi-Bellman case) or two parameters a, (3 
(in the case of Isaacs's equations). These parameters lie in index sets we will not 
display in the discussion. Thus we take as ingredients proper expressions of the 
form 



N ^9 N 

^U 



du 



(1.8) = - E <.(-)7w^r + E^"(-)^ + ^"W^(^) - /"( 



or 



dxidxj * dxi 



N ^2 N 

o u v-^ , du 



dxidxj ' ' dxi 

where all the coefficients are bounded with respect to a or a, j3. Hamilton-Jacobi- 
Bellman equations include those of the form 

(1.10) sup{£"u} = 

a 

while 



(1.11) 



supinf {r^-'^u} =0 

a 
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is a typical Isaacs's equation. The corresponding nonlinearities F have the form 
F{x,r,p,X) = sup [- tmce{A"{x)X) + {b"{x),p) + c"{x)r - /"{x)] 

a 

and 

F{x,r,p,X) =supinf[- tracc(A"'^^(a;)X) 

a P 

+ (6«''5(x),p>+c«''^(x)r-/«''5(x)], 

each of which is clearly also proper. Notice that in the first case F is convex in 
(r,p, X) while in the second case this is not so; indeed, if one allows for "unbounded 
envelopes" (i.e., coefficients that are unbounded in the parameters), one can show 
that essentially any proper F satisfying minor regularity assumptions can be rep- 
resented as a "sup inf of linear expressions as above. 

Indeed, the above process is quite general and does not require linear ingredients. 
Suppose F"''^ is proper for each a, /J. Then F{x, r,p, X) = sup„ inf/j F"'^{x, r,p, X) 
and F{x, r,p, X) = info, sup^ F"'^{x, r,p, X) are also evidently proper (for the mo- 
ment, wc set aside considerations of finitcness and continuity). 

Example 1.7. Obstacle and gradient constraint problems. A special case of the 
last remarks above is met in the consideration of "obstacle problems." Very general 
forms of such problems may be written 

(1.12) max{F{x, u, Du, D\), u - f{x)} = 
or 

(1.13) min{F(a;, u, Du, D'^u), u - f{x)} = 
or even 

(1.14) max{min{F(a;, u, Du, D'^u),u — f{x)}, u — g{x)} = 0. 

In accordance with remarks made in the previous example, if F is proper then so 
are (1.12)-(1.14). 

Likewise, "gradient constraints" may be imposed in this way. A typical example 
corresponds to max{F{x, u, Du, D'^u), \Du\ — g{x)} = 0. 

Example 1.8. Functions of the eigenvalues. For X G S{N) we let X\{X), . . . , \m{X) 

be its eigenvalues arranged in increasing order, Xi{X) <\i+i{X). li g{x,r,p, si, . . . , Sjv 
is defined on IR.^^+^ and is nondecreasing in r and each Sj, then F(x,r,p,X) = 
g{x,r,p,—Xi{X), . . . , —AAr(X)) is proper. For instance, F(X) = —max{Ai(X), .. . ,Aa 
'Xn{X), F{X) = -min{Ai(X),... ,Xn{X)} = -Xi{X) and F{X) = -(A2(X))3 
are degenerate elliptic. Another example is 

F{x,r,p,X) = -\ trace(X)|"'-^ trace(A:) + + c(x)r - f{x) 

where c > and m, g > 0. The corresponding equation is 

-|Aup-iAw + + c{x)u = f{x), 

which provides another example of the generality we are dealing with, even if we 
have no interpretation of this equation in mind. 
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Example 1.9. Sums and increasing functions of proper functions. If Fi is proper 
for i = 1, . . . , M, then so is Fi + • • • + Fm- More generally, if g{si, . . . , sm) is 
nondecreasing in each variable, then g{Fi, . . . , Fm) is proper. One may build very 
complex examples using the cases discussed above and these remarks. 

Example 1.10. Parabolic problems. We just observe that if {x, r, p, X) — > F{t, x, r, 
p,X) is proper for fixed t £ [0,T], then so is the associated "parabolic" problem 

(1.15) ut + F{t,x,u,Du,D'^u) =0 



when considered as an equation in the A'' + 1 independent variables {t,x). We 
mention only one example (there are, of course, infinitely many) that has some 
geometrical interest since it describes the evolution of a surface (given by a level set 
of the initial condition) with a motion along its normal with a speed proportional 
to the mean curvature 

(1.16) wt - div ( ) = 0. 



\Du\ 



Carrying out the differentiations yields 



^ ^ d^u du du 2 Q 

Ut u dxidxi dxi dxi ^ 

This may be written in the form (1.15) with 

F{x,p,X) = - trace (^(^I - X 



Example 1.11. Monge Ampere equations. The Monge-Ampere equation may be 
written as 

(1.18) u is convex, Aei{D^u) = f{x,u, Du) 

where f{x,r,p) > 0. We arc dealing here with the real Monge-Ampere equation, 
but everything that will be said adapts to the complex case and to other curvature 
equations. Allowing F to be discontinuous (even more, to become infinite), we may 
write (1.18) in our form by putting 



F{x,r,p,X) 



- det{X) + f{x,r,p) if X > 0, 
+00 otherwise; 



F is then degenerate elliptic. This follows from the fact that 



Hill Si if Si > 0, i = 1,... ,iV, 



—00 otherwise. 



is nondecreasing in each of its arguments and Example 1.8. 
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Example 1.12. Uniformly elliptic functions. This "example" is really a definition. 
One says that F{x,r,p,X) is uniformly elliptic if there are constants A, A > for 
which 

A tracc(P) < F(x,r,p, X - P) - F(x,r,p, X) < A trace(P) for P > 

and all x,r,p,X and then calls the constants A, A ellipticity constants. We note 
that L in Example 1.2 is uniformly elliptic with constants A, A exactly when A < 
Xi{A(x)) and Xj\[(A(x)) < A. One notes that sums of uniformly elliptic functions 
are again uniformly elliptic and that the sup inf process over a family of uniformly 
elliptic functions with common ellipticity constants produces another such function. 

Notes on §1. We will disappoint the reader in the following sections by not ap- 
plying the theory developed therein to the many examples given above. The goal of 
this section was to exhibit clearly the breadth and importance of the class of proper 
equations. We simply do not have enough space here to develop applications of the 
theory of these equations beyond that which follows immediately from the general 
results presented. 

Most of the examples listed have been considered via classical approaches. We 
give some references containing classical presentations: D. Gilbarg and N. S. Trudinger 
[81] is a basic source concerning linear and quasi-linear uniformly elliptic equations; 
O. A. Oleinik and E. V. Radkevic [139], J. J. Kohn and L. Nirenberg [106], and A. V. 
Ivanov [98] treat degenerate elliptic equations; W. H. Fleming and R. Rishcl [75], P. 
L. Lions [117], and N. V. Krylov [108, 109] are sources for Hamilton- Jacobi-Bellman 
equations; S. Benton [36] and P. L. Lions [116] discuss first-order Hamilton-Jacobi 
equations. Most of these references present some of the ways these equations arise. 

2. The notion of viscosity solutions 

It is always assumed that F satisfies (0.1) (i.e., F is proper) and, unless otherwise 
said, is continuous. To motivate the notions, we begin by supposing that u is 
(i.e., twice continuously differentiable) on and 

F{x,u{x),Du{x),D^u{x)) < 

holds for all x (that is, u is a classical subsolution of P = or, equivalently, a 
classical solution of P < in K^). Suppose that <^ is also and f is a local 
maximum of w — </j. Then calculus implies Du{x) = D(p{x) and D^u{x) < D^(p{x) 
and so, by degenerate ellipticity, 

(2.1) F{x, u{x), Dip{x), D'^<fi{x)) < F{x, u{x), Du{x), D^u{x)) < 0. 

The extremes of this inequality do not depend on the derivatives of u and so we 
may consider defining an arbitrary function u to be (some kind of generalized) 
subsolution of P = if 

(2.2) P(x, u{x), D(f{x), P'V(i)) < 

whenever (p is and i is a local maximum of u — (p. Before making any formal 
definitions, let us also note that u{x) < u{x) — ip{x) + ip{x) for x near x,ip € 
and Taylor approximation imply 

(2.3) u{x) < u{x) + (p,x — x) + ^{X{x — x),x — x) + o{\x — x\'^) as x ^ x 
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where p = Dip{x) and X = D'^ip{x). Moreover, if (2.3) holds for some {p,X) € 
X S{N) and u is twice differentiate at x, then p = Du{x) and D'^u{x) < X. 
Thus if u is a classical solution of F < it follows that F{x,u{x),p, X) < whenever 
(2.3) holds; we may also consider basing a definition of nondifferentiable solutions 
u of < on this fact. Roughly speaking, pursuing (2.2) leads to notions based 
upon test functions (p but does not immediately lead us, as will pursuing (2.3), to 
define ''{Du, D'^u)" for nondifferentiable functions u, which will turn out to be a 
good idea. For this reason, we begin by developing the line suggested by (2.3). Next 
we introduce a set O C on which F < is to hold and the appropriate notation 
to deal with inequalities like (2.3) "on O." At the moment, O is arbitrary; later we 
require it to be locally compact. Taking off from (2.3), if u: O ^ R, x G O, and 

(2.3) holds as O 5 X ^ x, we say {p,X) e Jq^u{x) (the second-order "superjet" 
of u x). This defines a mapping J^'^u from O to the subsets of x S{N). 

Example 2.1. By way of illustration, if u is defined on R by 

r for a; < 0, 

uix) = < , „ 

^ ' \ax+^x^ for a; > 0, 

then Jj^_:+ oju(O) = ((-oo, 0) x R) U ({0} x [0, oo)), while 

r ifa>0, 
j2.+y(o) = i {0} x [max{0,6},oo) if a = 0, 

[ ((a, 0) X M) U ({0} X [0, oo)) U ({a} x [b, oo)) if a < 0. 

Having thought through this example, the reader will see that J'q'^u{x) depends 
on O but realize it is the same for all sets O for which x is an interior point; 
we let J'^'^u{x) denote this common value. If we repeat the above discussion after 
switching the inequality sign in (2.3), we arrive at the definitions of the second-order 
"subjets" Jq~u, J'^'~u; equivalently, Jq~u{x) = —J^~^{—u){x), etc. 

We are ready to define the notions of viscosity subsolutions, supersolutions, and 
solutions. It will be useful to have the notations 

USC(O) = { upper semicontinuous functions u: O — > M}, 
LSC(O) = { lower semicontinuous functions u: O ^ M}. 

Definition 2.2. Let F satisfy (0.1) and O C M^. A viscosity suhsolution of F = 
(equivalently, a viscosity solution of F < 0) on O is a function u G USC(O) such 
that 

(2.4) F{x,u{x),p,X) <Q for all x&O and {p,X) & j'^'^u{x). 

Similarly, a viscosity supersolution of F = on O is a function u G LSC{0) such 
that 

(2.5) F{x,u{x),p,X)>Q for all xeO and {p,X) e j'^~u{x). 

Finally, « is a viscosity solution of F = in O if it is both a viscosity subsolution 
and a viscosity supersolution of F = in O. 
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Remarks 2.3. Since these "viscosity notions" are the primary ones for the current 
discussion, we immediately agree (at least, we hope you agree) to drop the term 
"viscosity" and hereafter simply refer to subsolutions, supersolutions, and solutions. 
This is a happy idea, as the term "viscosity," which lacks elegance, is an artifact of 
the origin of this theory in the study of first-order equations and the name was then 
motivated by the consistency of the notion with the method of "vanishing viscosity," 
which is irrelevant for many second-order equations. It follows from the discussion 
preceding the definition that, for example, if u is a solution of F < in C, is in 
a neighborhood of O, and u — ip has a local maximum (relative to O) at x G O, then 
(2.2) holds. Analogous remarks hold for supersolutions. These remarks motivate 
the requirement that a subsolution be upper semicontinuous, etc., in the sense that 
producing maxima of upper semicontinous functions is straightforward. Solutions, 
being both upper semicontinuous and lower semicontinuous, are continuous. One 
might ask if the validity of (2.2) for all (p E C"^ (with the maxima relative to O) for 
an upper semicontinuous function u is equivalent to u being a subsolution. This is 
so. In fact, if X € O then 

Jq~^u{x) = {{Dip{x), D'^ip{x)) : is and u — ip has a local maximum at f }; 

we leave the proof as an interesting exercise. 

We next record the definitions of the closures of the set- valued mappings needed 
in the next section. 

With the above notation, for x e O, we set 
(2.62 

/o'^u{x) = {(p, X) G X S{N) : 3(x„,p„, X„) G O x x S{N) 9 

{Pn,Xn) G j'^~^u{Xn) and {Xn,u{Xn),Pn,^n)-^{x,u{x),p,X)} 

and 
(2.7| 

Jo y-ix) = {{P,X) G X S{N) : 3(x-„,p„,X„) G O x x S{N) 5 

{Pn,Xn) G J^~u{xn) and {x„,u{xn),Pn,X„)^{x,u{x),p,X)}; 

we are abusing standard practice as regards defining closures of set- valued mappings 
a bit in that we put the extra condition u{x„) —>■ u{x) in the definitions while the 
graphs of the multifunctions j'^^u, u do not themselves record the values of 
u. The reader may note the use of expressions like " abbreviation 
for "the sequence a;„ satisfies lim„^ooa^n = x" etc. If a; G interior(O), we define 

2 + 2 — 

J ' u{x), J ' ^(a;) in the obvious way. 

Remark 2.4. If u is a solution of F < on O, then F{x,u{x),p, X) < for 

X G O and {p, X) G Jq^u{x). This remains true, for reasons of continuity, if 

2.+ 

{p,X) G Jq u{x) and F is continuous (or even lower semicontinuous). Similar 
remarks apply to supersolutions and solutions. 

Advice. We advise the reader to either skip the following material or to scan it 
lightly at the present time and proceed directly to the next section. The comments 
collected below can be referred to as needed. 

Remark 2.5. While the definitions above may seem reasonable, they contain sub- 
tleties. In particular, they do not define "operators" on domains in a familiar way. 
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For example, let 1 and consider the equation (m')^ — 1 = 0. The functions \x\ 
and — both satisfy this equation classically on R\{0}. However, the semijets of 
the function u{x) = — \x\ are given by 



so —\x\ is a solution of the equation, while if v{x) = \x\ then (0, 0) G ■P~v{Q) and 
V is not a solution. Similarly, one sees that docs solve — (w')^ + 1 = 0. Hence a 
solution oi F = Q need not be a solution of — F = (of course, if — F is not proper, 
as is typical when F depends on u,X, this is meaningless). It is also unknown if 
the information iJ, m, /, g £ C(R^) and u is a solution of H{Du) — f{x) = as well 
as H{Du) — g{x) = in implies f = g (although it is known that / = <? if either 
H is uniformly continuous or A'' = 1). 

We do not mention these things to disturb the reader; indeed, they will not 
explicitly appear in the theory. Our purpose is to emphasize that the theory to 
follow is not a variant of more classical developments. 

Remark 2.6. If we consider one of the simplest examples of a discontinuous upper 

scmicontinuous function, namely, the function given by u{x) = for x ^ 0, u{Q) = 
1, we have j2'+w(0) = R x R. The relation F(0, l,p,X) < for {p,X) € J^'+u{0) 
is then very restrictive. In fact, a primary role played by allowing semicontinuous 
subsohitions and supersolutons in the basic theory is only as a technical device to 
produce continuous solutions via Perron's method as described in §4. If solutions 
are produced by other methods (e.g., via numerical or other approximations) there 
may be no need to invoke results on semicontinuous functions. However, one of 
the simplest ways to show uniform convergence of approximation schemes, involves 
comparison of semicontinuous subsolutions and supersolutions. Thus one often 
deals with semicontinuity mainly as a device in proofs — as such, it is a powerful 
labor saving device. 

Remarks 2.7. We point out some facts concerning second-order semijets. 

(i) It is clear that J^~^u{x) is a convex subset of R-'^ x S{N); Example 2.1 shows 
that it is not necessarily closed. However, a little thought shows that if p € R''^, 
then {X : (p,X) € Jq^u{x)} is closed. 

(ii) If 9? e in some neighborhood of O, then Jq~^{u—ip){x) = {{p—Dip{x),X— 
D'^ip{x)) : {p,X) G J'q'^u{x)}. As a consequence, the same statement holds if Jq'^ 
is replaced everywhere by Jq~ , Jq^ , or Jq . Indeed, if {q,Y) e Jq^{u — </?)(x). 



u{x) —ifi{x) < u{x) — ip{x) + {q,x — x) + ^{Y{x — x), X — x) + o{\x — xl"^) as x ^ x 



' {(-1,X) : X > 0} if a; > 0, 

j2'+u(a;)=i ((-1,1) X R) U ({-1,1} X [0,oo)) ifa; = 0, 
^ {(1,X) : X > 0} if a; < 



and 




then 



(2.8) 



and 



^p{x) — ip{x) = {Dip{x), {x — x)) + ^{D'^ip{x){x — x),x — x) + o{\x — x^) 
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imply that {q + Dip{x),Y + D'^(p{x)) G Jq~^u{x) and so Jq'^{u — (p){x) C {{p — 
Dip{x),X — D^ip{x)) : {p,X) G Jq^u{x)}. The other inclusion follows from this as 
well, since Jq^u{x) = J'q^{{u — </?) + ip){x). It is also clear that one always has 

Jq^{u + v){x) D Jq~^u{x) + Jq~^v{x). 

(iii) We consider jQ~^ifi{x) when ip G C'^{M.^) and will end up with a general 
statement corresponding to Example 2.1. In view of (ii), we may as well assume 
that ^ = 0, and we will write "Zero" for the zero function. We know that if 
X e interior(e)), then J^~^ZeTo{x) = J'^'+Zevo{x) = {(0,X)} : X > 0}. In general, 
{p,X) e J^'+Zero(^) if 

(2.9) Q<{p,x-x) + \{X{x-x),x-x)-\-o{\x-x\^) bsO^x^x. 

Assuming that x„ € O, < | I — » 0, and {xn — x)/\xn — x\ q, we may put 

X = a;„ in (2.9), divide the result by | I and pass to the limit to conclude that 

(2.10) < (p, q) for q G UTo(i) 

where UTei(i;) is the set of "generalized unit tangents" to O at x; it is given by 

(2.11) UTo(f) = \q : 3a;„ e 0\{x}, Xn ^ x, and -> 

\^ \Xn X\ 

If is a smooth A^-submanifold of with boundary and x G dO, then the 
generalized tangent cone 

To{x)= convex hull(UTo(i)) 

is a halfspacc and O has an exterior normal n at x. In this event, (2.10) says that 
p = —Xn for some A > 0. Moreover, if p = 0, (2.9) then implies that < {Xq, q) 
for q e VToix), and we conclude that 

(2.12) (0, X) € Jo'+Zero(x) if and only if < X 

provided To{x) is a halfspace. 

Life is more complex if p = —An and A > 0. In this case, (2.9) reads 

(2.13) X{n,x — x) — ^{X{x — x),x — x) < o{\x — x\'^) as C 9 a; — > x. 
We study (2.13) when x = and we can represent O near in the form 

(2.14) {{x,xn) - Xn <9{x)} 

where x = {xi,. . . ,xjv-i), g{x) = \{Zx,x) + o(|xp)}, and Z G S{N — 1). That is, 
we assume that the boundary of O is twice differentiable at and rotate so that the 
normal is n = e^r = (0, ... ,0, 1). With these normalizations, we put {x,g{x)) into 
(J.13) to find A(Zi,i) - (X(i,0),(i,0)) < o(|i|2) or \Z < P^-iXPn-i where 
Z{x, Xn) = {Zx, 0) and Pjv-i is the projection on the first N —1 coordinates. It is 
not hard to see that this is also sufficient. 
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When unraveled, the above considerations lead to the following conclusion. Let 
O be an A^-submanifold of with boundary, x e 90, dO be twice differentiable 
at X, ft be the outward normal, Tqo{x) be the tangent plane to dO at x, and 
P G S{N) be orthogonal projection on Tqo{x) (which we regard as a subspace 
of R^). Finally, let S be the symmetric operator in Tao{x) corresponding to the 
second fundamental form of dO at x (oriented with the exterior normal to O) 
extended to IR^ by 5n = 0. Then 
(2.15) 

f (p, X) G J(p'^(^(a;) if and only if either p = Dip{x) and D'^(p{x) < X, or 
\p = Difix) - An, A > and PD'^ip{x)P < PXP - XS. 

Noting that Pn®nP = 0, we see that if {p,X) e Jq'^u{x), 5* < 0, and A > 0, then 
(p — An, X + ij,n®n) € Jq~^u{x) for fi G M.. For later use, we denote by Xx the 
particular choice of X such that {Dip{x) — An, X) G jQ'^ip[x) given by 

(2.16) XA(C + /3n) = r'V(*)(^ + /3n) + A5C for ^ G Tao(£),/3 e M. 
Before leaving this topic, we note that the set of exterior "normals" 

(2.17) N{x) = (n G : (n, g) < for g G Toix)) 
can also be described by 

(2.18) N{x) = {n&m.^ : {n,x-x) < o{\x - x\) bs O 3 x ^ x) 
no matter how regular dO is at x and that we showed above that 

(2.19) Jo'^vix) C {Dif{x) - N{x)) X S{N). 

Notes on §2. Viscosity solutions were introduced by M. G. Crandall and P. L. 
Lions [51] (following an announcement in [50]); analogies with S. N. Kruzkov's 
theory of scalar conservation laws [107] provided guidance for the notion and its 
presentation. The presentation of [51] emphasized the first-order case, owing to 
the fact that uniqueness results were only obtainable for this case at the time. It 
was also pointed out in [51] that there were several equivalent ways to formulate 
the notion of viscosity solutions; one of these was intimately connected with earlier 
work by L. C. Evans [64, 65] concerning "weak passages to the limit" in equations 
satisfying the maximum principle. In these developments, directions were indicated 
by aspects of nonlinear functional analysis and nonlinear semigroup theory. 

It became apparent that working with these alternative formulations was superior 
to the approach taken in [50], and perhaps M. G. Crandall, L. C. Evans, and P. 
L. Lions [47] is the first readable account of the early theory. (We note that the 
definition most emphasized in [47] is not the one put foremost here.) 

The necessity of defining the semijets on closed sets became apparent with the 
investigation of boundary value problems. We return to this issue in §7, but a couple 
of significant papers here are P. L. Lions [121] and H. M. Soner [149] (see also M. 

G. Crandall and R. Newcomb [58]). The closure of the semijets was introduced in 

H. Ishii and P. L. Lions [96] and was fully exploited in M. G. Crandall [50] and M. 
G. Crandall and H. Ishii [48]. Analogous operations have long been employed by 
"nonsmooth analysts." 
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The fact that the theory required only appropriate scmicontinuity was realized 
earlier on, but the first striking applications of the fact (to existence, stability, 
optimal control, . . . ) were given in H. Ishii [86, 90], G. Barles and B. Perthame [25, 
26]. 

Let us finally mention that some structural properties of viscosity solutions were 
and will be ommitted from our presentation and, in addition to the above references, 
some can be found in P. L. Lions [124], R. Jensen and P. E. Souganidis [105], and 
H. Frankowska [78]. 

3. The maximum principle for semicontinuous functions 
and comparison for the dirichlet problem 

Let f2 be a bounded open subset of M^, f2 be its closure, and dQ be its boundary. 
Suppose u e USC(f2) is a solution of F < in fl, v G LSC(f2) is a solution of F > 0, 
and u < V on dfi. We seek to show that u < v on Q. 

In the event u and v are classical sub and supersolutions, we could employ the 
classical "maximum principle." Let us recall some elementary facts in this regard. 
Let w be defined in a neighborhood of f e R^. If there exists (p, X) G x S{N) 
such that 

(3.1) w{x) = w{x) + {p,x — x) + ^{X{x — x),x — x) + o{\x — x]"^) as a; — > f , 

we say that w is twice differentiable at x and Dw{x) = p, D'^w{x) = X (it being 
obvious that p and X are unique if they exist). It is clear that w is twice differ- 
entiable at X if and only if J^iv{x) = J^'+w(x) n J^'^w(x) is nonempty (in which 
case J^w(x) = {{Dw{x), D'^w{x))}). Classical implementations of the maximum 
principle are based on the fact, already used above, that if w is twice differentiable 
at a local maximum x, then Dw{x) = and D^w{x) < 0. Thus, if u and v are 
twice clifFcrentiable everywhere and w = u — v has a local maximum x E fl, we 
would have Du{x) = Dv{x) and D'^u{x) < D'^v{x) and then, in view of (0.1), 

F{x,u{x),Du{x),D'^u{x)) < < F{x,v{x), Dv{x), D'^v{x)) 

< F{x, v{x), Du{x), D^u{x)). 

In the event that F{x, r, p, X) is strictly nondecreasing in r (a simple but illustrative 
case) , it follows that u ~ v is nonpositive at an interior maximum and so u < in 
n since u < v holds on dfl. 

We seek to extend this argument to the case u G USC(f7), v G LSC(il). We 
are unable to simply plug {Du{x), D'^il{x)) and {Dv{x), D'^v{x)) into F since these 
expressions must be replaced by the set-valued functions J^'^u and J^'^u (and 
their values may well be empty at many points, including maximum points of 
u{x) — v{x)). To use J'^'~^u and J'^'~v, we employ a device that doubles the number 
of variables and then penalizes this doubling. More precisely, we maximize the 
function u{x) — v{y) — (a/2)|a; — over fi x fi; here a > is a parameter. As 
a — > oo, we closely approximate maximizing u{x) — v{x) over More precisely, 
we have 

Lemma 3.1. Let O be a subset ofR^, u G USC(O), v G LSC(C>) and 



(3.2) 



M„ = sup {u{x) - v{y) - fix - y\^) 

OxO 
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for a > 0. Let Ma < oo for large a and {xa, Va) be such that 



(3.3) 



a 



lim {Ma - iu{Xa) - v{ya) - f - ^aP)) = 0. 



Then the following holds: 



(3.4) 



< 



(i) liniQ^oo "Ixq - 2/a|^ = and 

(ii) lim^j^oo Ma = ^{x) — v{x) = suPq{u{x) — v{x)) 
whenever x € O is a limit point of Xa as a ^ oo. 



Deferring the elementary proof of the lemma and returning to our sub and su- 
persolutions u G USC(f2), v G LSC(O) satisfying u < v on dfl, we note that 
Ma = snpQ^Q{u{x) — v{y) — {a/2)\x — yp) is finite since ^(a;) — v{y) is upper 
scmicontimious and is compact. Since wc seek to prove u < v, we assume to the 
contrary that u{z) > v{z) for some z G fi; it follows that 



Choosing (xa,y„) so that M„ = u{xa) — v{ya) — {a/2)\xa — Val^ (the maximum is 
achieved in view of upper semicontimiity and compactness), it follows from (3.4) (i), 
(ii) and u < v on dfl that {xa,ya) G x for a large. The next step is to use the 
equations to estimate M„ and contradict (3.5) for large a. This requires producing 
suitable values of J'^''^u and J^'~?;, and wc turn to this question. 

To know what to look for, let us proceed more generally and assume that u, v are 
defined in neighborhoods oi x,y G and twice differentiable at ;c, y respectively. 
Assume, moreover, that cp is near (x, if) in K^xM^ and (i, if) is a local maximum 
of u{x) — v{y) — (p{x, y). Applying the classical maximum principle to this situation 
(in the 2N variables (a;, y)), we learn that Du{x) = Dx(p{x, y), Dv{y) = —Dy(p{x, y), 
and 



where X = D^u{x), Y = D'^v{y). Note that with the choice ip{x, y) = {a/2)\x-y\'^ 
the above reads 



where / will stand for the identity matrix in any dimension and, since the right- 
hand side annihilates vectors of the form ( ^ ) (also written *(^, ^) where *Z denotes 
the transpose of a matrix Z), (3.7) implies X < Y, making further contact with 
the maximum principle. 

It is a remarkable fact that perturbations of the above results may be obtained 
in the class of semicontinuous functions. The main result we use in this direction 
is the following theorem, which is formulated in a useful but distracting generality. 
(See (2.6) and (2.7) regarding notation below.) 



(3.5) 



Ma > u{z) - v{z) = 6>0 for a > 0. 



(3.6) 




(3.7) 
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Theorem 3.2. Let Oi be a locally compact subset of R'^' for i = 1,. . . ,k, 

Ui G USC(Oi), and (p be twice continuously differentiable in a neighborhood of O. 
Set 

w{x) = ui{xi) H \-Uk{xk) for X = {xi, ■■■ ,Xk) 

and suppose x = {xi, . . . , Xk) G O is a local niMximum. ofw — (p relative to O. Then 
for each e > there exists e >S(iVj) such that 

(DxMx), Xi) e Jo. Ui{xi) fori = l,... ,k, 
and the block diagonal matrix with entries Xi satisfies 

'Xi ... 

(3.8) -{i + ||A||V< ( ; ••. M <A + eA' 



£ 



... Xk. 



where A = D'^ip{x) € S{N), N = Ni + --- + Nk. 

The norm of the symmetric matrix A used in (3.8) is 

ll^ll = sup{|A| : A is an eigenvalue of A} = sup{\{A^,0\ ■ ICI < !}• 

We caution the reader that this result is an efficient summarization of the analytical 
heart of the theory suitable for the presentation we have chosen and its proof, which 
is outHned in the appendix, is deeper and more difficult than its applications that 
we give in the main text. See also the notes to this section. 

In order to apply Theorem 3.2 in the above situation, we put k = 2, Oi = O2 = 

2 2 -|- 

il, ui = u, U2 = —V, ip{x,y) = {a/2)\x — and recall that Jq v = —Jq {—v). 
In this case 



Dx<f{x,y) = -Dyf{x,y) = a{x-y), A = a(^Jj ^ 
A^ = 2aA, \\A\\ = 2a, 

and we conclude that for every s > there exists X, Y G S{N) such that 

(3.9) ia{x - y), X) e Ja'^u{x), {a{x - y), Y) e J^'viy) 

and 

'I . ^ \ f I 0\ ^ f X \ -/ 
e 



2") 1 /) ^ ( _y <"(l + 2ea), _^ J 



Choosing s = 1/a yields the elegant relations 
(3^10) -Sa(' »1<(-J »)<3a^ ' 



I - \ -Y - \ -I I 
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in all, wc conclude that at a maximum (x, y) of u{x) — v{y) — {Q/2)\x—y\^ there exists 
X, — Y G S{N) such that (3.9) and (3.10) hold. Observe that, in comparison with 
(3.7), we have worsened the upper bound by a factor of 3 (but can still conclude from 

(3.10) that X <Y) and obtained a lower bound. We will not actually use the lower 
bound in what follows, but its presence corresponds to an essential compactness 
issue in the proof of the theorem. 

Let us finish this chapter of the story. Writing (f,y) in place of {Xa,ya) for 
simplicity, we only need to assume that 

(3.11) F{x, u{x), a{x-y),X)<0< F{y, v{y), a{x -y),Y) 

is inconsistent with the other information we have in hand, which includes 



(3.12) 



< ^ < Ma = u{x) - viy) - {a/2)\x - < u{x) - v{y) and 
a\x — — > as a — > 00 



(the second assertion follows from (3.4)(i)) and (3.10). For example, assume that 
there exists 7 > such that 

(3.13) 7(r - s) < F{x,r,p,X) - F{x,s,p,X) for r > s, {x,p,X) eTlxR^ x S{N) 

and there is a function w : [0, 00] [0, 00] that satisfies a;(0+) = such that 

(3 14") / ^^y''^''^^^~y^'^^ ~ F{x,r,a{x -y),X) < u;{a\x - y\^ + \x - y\) 
\ whenever x,y£n, r e R, X, F e S{N), and (3.10) holds. 

Proceding, we deduce from (3.12), (3.13), (3.11), and (3.14) that 

l5 < l{u{x) - v{y)) < F{x, u{x), a{x - y), X) - F{x, v{y), a{x - y), X) 

= (F{x,u{x),a{x -y),X) - F{y,v{y),a{x -y),Y)) 
(3-15) + (^(^^ ^(^)^ -y),Y)- F{x, v{y), a{x - y), X)) 

<uj{a\x-y\^ + \x-y\)\ 

here wc used (3.11) to estimate the first term on the right by and (3.14) on the 
second term. Since uj{a\x — + |x — y|) — >■ as a — > 00 by (3.12), we have a 
contradiction. We have proved 

Theorem 3.3. Let Q be a bounded open subset ofR^, F e C(f2 x M x x S{N)) 

be proper and satisfy (3.13), (3.14). Let u £ USC(ri) [respectively, v S LSC(ri)) be 
a subsolution {respectively, supersolution) of F = in fl and u < v on dfl. Then 
u < V in Q. 

Remark 3.4. We motivated the structure condition (3.14) by the natural way it 

arose in the course of analysis and (3.14) was used in the proof and not degenerate 
ellipticity. However, (3.14) implies degenerate ellipticity. To see this, suppose that 

(3.14) holds, X,Y G S{N) and X <Y. Observe that then for ^, r? G M^, and e > 

{X^, - {Yv, V) < (n, - {YV, V) - 2(rr?, ^ - r?) + (F(? - r;), ^ - 7?) 
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The last inequality may be written 

Thus if a > (l/3)max{||X|l, ||F||, (1 + ||y||/e)||F||} and e is sufficiently small, the 
pair X, y + e/ satisfy (3.10). Next we fix x, r, and p, and put y = x—p/a in (3.14) 
so that by (3.14) 

F(x-^,r,p,Y + el)- F{x, r,p, X) < (^^{\p\^ + \p\)j 

and then we let a — > oo and £ | and conclude that F is degenerate elliptic. 

Remark 3.5. The proof of Theorem 3.3 adapts to provide modulus of continuity 
estimates on solutions of F = 0. 

Example 3.6. We turn to some examples in which (3.14) is satisfied. First 
notice that (3.14) evidently holds (with u the modulus of continuity for /) for 
F{x, r,p, X) = G{r,p, X) — f{x) if G is degenerate elliptic and / is continuous; this 
is because (3.10) implies X <Y. 

Secondly, the linear first order expression {b{x),p) satisfies (3.14) if 

a{b{y) - b{x),x - y) < ui{a\x - yf + \x- y\) 

for some ut and this holds with oj{r) = cr if there is a constant c > such that 
{b{x) — b{y),x — y)> —c\x — t/p, i.e., the vector field b + cl is "monotone." In fact, 
it is not hard to see that this is a necessary condition as well. 
Next, the linear expression 

(3.16) G{x, X) = - trace(*S(x)S(x)X), 

where S maps Q into the N x N real matrices, is degenerate elliptic (it is a special 
case of Example 1.2) with A{x) = *I](x)S(x)). We seek to estimate G{y,Y) — 
G{x,X) when (3.10) holds. Multiplying the rightmost inequality in (3.10) by the 
nonnegative symmetric matrix 

/*S(x)E(^) *E(2/)S(x)) 



and taking traces preserves the inequality and yields 

(3 17) " " trace(*E(a;)S(.T)X - 'j:{y)E{y)Y) 

< 3a trace((*S(x) - *S(y))(S(a;) - S(y)), 

so if E is Lipschitz continuous with constant L then 

G{y, Y) - G{x, X) < 3L'^a\x - yf 

and we may choose cj(r) = 3i^r. We note that if r(a;) : O — > S{N) is Lipschitz 
continuous and r(a;) > in Q, then 

- trace(r(a;)X) = - trace(r(a;)^/^r(a;)^/^X) 
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is an example of the sort just discussed since T^xY^"^ is also Lipschitz continu- 
ous. Finally, it is known that it r(a;) > and T G W'^'°°, then T^/^ is Lipschitz 
continuous. 

Now we note that the processes of forming sums or sups or sup infs, etc., above 
produce examples obeying (3.14) if the ingredients obey (3.14) with a common oj. 
Finally, we may produce examples satisfying (3.13) and (3.14) by putting 

F{x,r,p,X)='yr + G{x, r, p, X) 

whenever G is degenerate elliptic and satisfies (3.14). 

We deferred discussion of the elementary Lemma 3.1 since it is better presented 
in a generality not required above. Indeed, Lemma 3.1 may be obtained from the 
following proposition via the correspondences 2A'' — > M, O x O ^ O, (a;, y) — > 
X, u{x) - v{y) ^{x), {l/2)\x - ^ *(ar). 

Proposition 3.7. Let O be a subset ofR^, $ e USC(O), * G LSC(0), * > 0, 
and 

(3.18) M„ = sup($(a;) - a^{x)) 

o 

for Of > 0. Let — oo < lima^cc < oo and Xa ^ O be chosen so that 

(3.19) lim {Ma - mxa) - a*(x„)) = 0. 

a— ^oo 

Then the following hold: 

(i) lima-^oo a'^{xa) = 0, 



(3.20) 
Proof. Put 



(ii) *(x) = and lima^oo Ma = $(i) = sup{^(2,)^o} ^i^) 

whenever x & O is a limit point of Xa as a ^ oo. 

6a=Ma- ($(Xa) - a'^iXa)) 



SO that 5a as a — » 00. Since ^' > 0, Ma decreases as a increases and 
limc^oo Ma exists; it is finite by assumption. Moreover, 

Ma/2 > HXa) - -^-^{Xa) > ^(a^a) " a^{Xa) + '^'^{Xa) > - (5„ + -*(Xa), 

so 2{Ma/2 — Ma + (5a) > a'^{xa), which shows that a'i'{xa) — * as a ^ oo. 
Suppose now that a„ ^ oo and Xa^ — > x G C Then "^{xa^) and by the lower 
semicontinuity ^'(i) = 0. Since 

^{Xan) - a„*(Xa„) > Ma„ - 5a„ > SUp <^{x) - 5a„ 

{*=0} 

and $ is upper semicontinuous, (3.20) holds. 

Remark 3.8. We record, for later use, some observations concerning maximum 
points (:r, y) of u{x) — v{y) — ip{x — y) over O x for other choices of ip besides 
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{a/2)\x\'^ and the imphcations of Theorem 3.2. We are assuming that u is upper 
semicontinuous, v is lower semicontinuous, and if G C"^. In this application, the 
matrix A of Theorem 3.2 has the form 

(3.21) ^={-z ~z)' "^^^^^ ^ = - y), 



so 



(3.22) = 



2Z2 -2Z2 
-2^2 2^2 



and \\A\\ =2||Z||. 

Choosing e = in (3.8), we conclude that there are X, Y £ S{N) such that 

(3.23) {D^{x -y),X)G JI+u{x), {Dip{x -y),Y)€ J^'viy), 

and 



(3 24) -2m (' ^Ur^ Z + j^Z^ -iZ+Jikl^' 

(3.24) 2\\A\\^^ ^)<^^ _Y ) ^ [-iZ + ^Z-) Z + ^Z^ 



in particular, 

(3.25) ||X||, ||F|| < 2||A|| and X <Y. 

Notes on §3. This section is self-contained except for the proof of Theorem 3.2, 
which is due to M. G. Crandall and H. Ishii [48]. This result distills and sharpens 

the csscince of a Une of development that runs through several of the references 
listed below; the proof is explained in the appendix for the reader's convenience. 
We want to emphasize at the outset that our presentation of the proof of the 
comparison theorem, based as it is on Theorem 3.2, is but one among several 
possibilities. Moreover, other approaches may be useful in other situations. For 
example, one may use elements of the proof of Theorem 3.2 directly and look at 
corresponding manipulations on solutions — this approach seems to be the; leading 
presentation for the infinite-dimensional extensions of the theory and allows for 
adaptations on the test functions and the notion itself and avoids the "semijets" 
or generalized derivatives. Another presentation selects another building block, 
namely, regularization by supconvolutions (which occur in the proof of Theorem 
3.2) and its effect on viscosity solutions; this procedure is especially helpful in some 
particular problems (integrodifferential equations, regularity issues) and stresses 
the regularization procedure. This procedure has some analogues with the use of 
mollification in the study of linear partial differential equations. With apologies for 
confusing the reader, we are saying that a full grasp of all the elements of the proof 
of Theorem 3.2 and exposure to other presentations (e.g., R. Jensen [102] and H. 
Ishii and P. L. Lions [96]) may prove valuable. 

The first uniqueness proofs for viscosity solutions were given for first-order equa- 
tions in M. G. Crandall and P. L. Lions [51] and then M. G. Crandall, L. C. Evans, 
and P. L. Lions [47] . The second-order case remained open for quite a while during 
which the only evidence that a general theory could be developed was in results for 
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Hamilton-Jacobi-Bellman equations obtained by P. L. Lions [119, 118]. The proof 
in these works involved ad hoc stochastic control verification arguments. 

A breakthrough was achieved in the second-order theory by R. Jensen [101] 
with the introduction of several key arguments; some of these were simplified in R. 
Jensen, P. L. Lions, and P. E. Souganidis [104], P. L. Lions and P. E. Souganidis 
[132], and R. Jensen [102]. In particular, the use of the "supconvolution" regu- 
larization (see the appendix), which is a standard tool in convex and nonsmooth 
analysis, was somewhat "remise au gout dujour" by J. M. Lasry and P. L. Lions 



Progress in understanding these proofs so as to be able to handle more examples 
was made by H. Ishii [89] who introduced matrix inequalities of the general form 
(3.7). This work contains an example (following [89, Theorem 3.3]) showing the 
optimality of condition (3.14) in the sense that we cannot replace the right-hand side 
of (3.14) by, for example, uj{a\x — yf + \x — y\) with 9 <2 (but see §5. A). Estimating 
the left-hand side of (3.14) differently, for example in the form ag{\x — y\) + \x — y\, it 
is possible to prove uniqueness in cases in which g{r)/r'^ is unbounded near r = by 
refining the arguments of [49]. A more complete understanding together with the 
sharpest structure conditions (up to the present) were achieved in [96], which also 
contains many examples. We mention that the structure condition (3.14) exhibited 
in this section pays no attention to whether the equation is of first order or possess 
stronger ellipticity properties. In the case of uniformly cUipitic F, much more can 
be done [96]. See also Jensen [102] as regards structure conditions. 

A useful sharpening and improved organization of the analytical essence of the 
theory were contributed by [46]. The presentation here was based on Theorem 3.2; 
as remarked above, it follows from [48] and the generality and proofs presented 
there represent the current state of the art (see the appendix). 



Let be an arbitrary open subset of M^. By a solution (respectively, subsolution, 
etc.) of the Dirichlet problem 



we mean a function u G C{fl) (respectively, u e USC(n), etc.) that is a (viscosity) 

solution (respectively, subsolution, etc.) of F = in f7 and satisfies u{x) = 
(respectively, u{x) < 0, etc.) for x G dft. (We note that this formulation imposes 
the boundary condition in a strict sense that will be relaxed in §7.) 

Recall that wc always assume that F is proper and, unless otherwise said, con- 
tinuous . To discuss Perron's method, we will use the following notations: if 
u: [-00, do] where O C M^, then 



[110]. 



4. Perron's method and existence 



(DP) 




'u) =0 in f2, 
u = on dfl 



(4.1) 



u*{x) = limrj^o sup{u(y) : y G O and |y — a;| < r}, 
u*(x) = lim,.|o inf {7i(y) : y G O and \y — x\ < r}. 



One calls u* the upper semicontinuous envelope of m; it is the smallest upper semi- 
continuous function (with values in [—00,00]) satisfying u < u* . Similarly, is 
the lower semicontinuous envelope of u. 
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Theorem 4.1. (Perron's Method). Let comparison hold for (DP); i.e., if w is a 
suhsolution of (DP) and v is a supersolution of (DP), then w < v. Suppose also 
that there is a suhsolution u and a supersolution u of (DP) that satisfy the boundary 
condition u^{x) = u*{x) = for x G dO,. Then 

(4.2) W{x) = sup{w{x) •.u<w<u and w is a suhsolution of (DP)} 

is a solution of (DP). 

Theorem 3.3 provides conditions under which comparison holds for (DP). The 
proof consists of two steps. The first one is 

Lemma 4.2. Let O C &e locally compact, F e LSC(C' x K x x S{N)), 
and J- be a family of solutions of F < in O. Let w{x) = sup{u{x) : u € J-} and 
assume that w*{x) < oo for x € O. Then w* is a solution of F <0 in O. 

Allowing F to be merely lower semicontinuous does not affect the proof and is 
used later. The information produced while proving this lemma is needed later as 
well, so we isolate the essential point in another result. 

Proposition 4.3. Let O C be locally compact, v G USC(O), z e O, and 
{p,X) G j'^^v{z). Suppose also that Un is a sequence of upper semicontinuous 
functions on O such that 



(4.3) 

Then 

(4.4) 



(i) there exists x„ G O such that (.t„, w„(a;„)) ^ (z,w(z)), 

(ii) if Zn € O and Zn ^ x € O, then limsup„^(^ Un{zn) < v{x). 



there exists f„ G O, (p„,X„) G j'^^Un{xn) 

such that iXn,Un{Xn),Pn,Xn) (z , v{z) , p, X) . 



Proof. Without loss of generality we put z = 0. By the assumptions, for every 
5 > there is an r > such that = {x G O : \x\ < r} is compact and 

(4.5) v{x) <v{0) + {p,x) + ^{Xx,x) +S\xf for x e Nr. 

By the assumption (4.3)(i), there exists a;„ G O such that (a;„, M„(a;„)) {0,v{Oj). 
Let Xn G iVr be a maximum point of the function Un{x) — {{p,x) + ^{Xx,x) + 251x1"^) 
over Nr so that 

Un{x) < Un{Xn) + {p,X-Xn) + h{{Xx,x) - (Xf „, f „)) + 25(|a;|^ - \Xnf) 

(4.6) 

for X G Nr. 

Suppose that (passing to a subsequence if necessary) .i„ ^ y as n — > oo. Putting 
a; = x„ in (4.6) and taking the limit inferior as n — > oo, we find 

v{0) < liminf u„(x„) - (p,y) - \{Xy,y) - 25\y\^\ 

on the other hand, by (4.3)(ii) Hminfu„(x„) < v{y) while (4.5) implies v{y) — 
{p,y) - \{Xy,y) - 25\y\'^ < v{0) - 6\y\^. We conclude that 

v{0) < liminf w„(i„) - {p,y) - \{Xy,y) - 25\y\'^ 

< viy) - - l{Xy,y) - 2S\yf < v{0)-S\yf. 
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Prom the extreme inequalities we learn y = 0, so i„ ^ (without passing to 
a subsequence), and then from the first inequality and (4.3) (ii) one sees that 
v(0) = lim„^oo Mn(^n)- Since we have {p + 4(5i„ + Xxn,X + ASI) e jQ^Un{xn) 
for large n, we are nearly done. To conclude, we merely note that the set of 
{q,Y) e X S{N) such that there exists z„ e O, (p„,X„) G J^^Unizn) such 
that {zn,Un{zn),Pn, Xn) {0,v{0), q,Y) IS closcd and contains {p,X + 461) for 
5 > by the above. 

Proof of Lemma 4.2. With the notation of the lemma, suppose that z E O and 
{p,X) e Jq~^w*{z). We seek to show that F{z,w*{z),p,X) < 0. It is clear that we 
may choose a sequence (x„, u„) G O x such that (x„, u„(x„)) — > (z, w*{z)) and 
that (4.3) then holds with v = w* . Hence, by the existence of data satisfying (4.4) 
and the fact that each u„ is a subsolution, we may pass to the limit in the relation 
F{xn,Un{xn),Pn,Xn) < to find F {z , w* {z) , p, X) < as desired. 

The second step in the proof of Theorem 4.1 is a simple construction that we 
now describe. Suppose that f2 is open, u is a solution of F < 0, and u* is not a 
solution of F > 0; in particular, assume G O, and we have 

(4.7) F(0,u*(0),p,X) < for some (p,X) e J^"w*(0). 

Then, by continuity, us^^{x) — -^,(0) + S + {p,x) + ^{Xx,x) — 7|a;p is a classical 
solution of < in i?r = {a; : \x\ < r} for all small r, (5, 7 > 0. Since 

u{x) > u»{x) > Ut:{0) + {p,x) + ^{Xx,x) + o{\x\'^), 

if we choose 6 = (r^/8)7 then u{x) > us,^{x) for r/2 < \x\ < r if r is sufficiently 
small and then, by Lemma 4.2, the function 

[/•(j.) _ / max{u(a;),U5,^(a;)} if|a;| < r, 
\ u{x) otherwise, 

is a solution of < in O. The last observation is that in every neighborhood of 
there are points such that U{x) > u{x); indeed, by definition, there is a sequence 
{xn,u{xn)) convergent to (0,u*(0)) and then 

lim {U{xn) - u{xn)) = "5,7(0) - u*(0) = u*(0) + "^O) > 0- 

n— »oo 

We summarize what this "bump" construction provides in the following lemma, the 
proof of which consists only of choosing r, 7 sufficiently small. 

Lemma 4.4. Let Q be open and u he solution of F < in fl. If fails to be 
a supersolution at some point x, i.e., there exists {p,X) G u^{x) for which 
F{x,u^{x),p,X) < 0, then for any small k> there is a subsolution of F <0 
in Cl satisfying 



(4.8 



Uk{x) > u{x) and supQ(f7K — m) > 0, 
Uk{x) = u{x) for x G Q,,\x — x\ > k. 



Proof of Theorem 4.1. With the notation of the theorem observe that < W* < 
W <W* <u* and, in particular, W^, = W = W* = on dQ. By Lemma 4.2 
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W* is a subsolution of (DP) and hence, by comparison, W* < u. It then follows 
from the definition of W that W = W* (so W is a subsolution). If W* fails to be 
a supersolution at some point x G fi, let Wi^ be provided by Lemma 4.4. Clearly 
u < Wk and = on dil for sufficiently small k. By comparison, W^, < u and 
since W is the maximal subsolution between u and u, we arrive at the contradiction 
< W. Hence W* is a supersolution of (DP) and then, by comparison for (DP), 
W* = W < W* , showing that W is continuous and is a solution. 

Remarks 4.5. The assumption that f7 was open in Lemma 4.4 was used only to 
know that classical subsolutions in small relative neighborhoods of points of were 
subsolutions in our sense. In order to generalize this and to formulate the version 
of Theorem 4.1 we will need in §7 (which we did not do above for pedagogical 
reasons), we now make some remarks and invite the reader to ignore them until 
they are called for later. Suppose O is locally compact, G+, G_ are defined on Ox 
R X X S{N) and have the following properties: G'^ is upper scmicontinous, G- 
is lower semicontinuous, and classical solutions (twice continuously differentiable 
solutions in the pointwise sense) of G+ < on relatively open subsets of O are 
solutions of G_ < 0. Suppose, moreover, that whenever u is a solution of G- < 
on O and v is a solution of G+ > on O we have u <v on O. Then we conclude 
that the existence of such a subsolution and supersolution guarantees that there is 
a unique function u, obtained by the Perron construction, that is a solution of both 
G-i- > and G_ < on C The proof is unchanged except in trivial ways. 

Example 4.6. Theorem 4.1 leaves open the question of when a subsolution u and 
a supersolution u of (DP) that vanish on dO. can be found. Let us consider this 
problem for the equation 



where G is degenerate elliptic. Sometimes there may be an obvious choice for u or 
m; e.g., if 



then w = is a subsolution. We assume (4.10) and seek a supersolution "near 90' 
in the form 



(4.9) 




'u) = in ri, 
u = on do. 



(4.10) 



G(.T,0,0)<0 for a; e 17, 



(4.11) 



tti(x) = M(l - e-^'^^^)) 



where M, A > are parameters to be chosen later and 



d{x) = inf \y - x\ 

yedQ 



is the distance to dCl. For 7 > set 




if dO, is of class G^ and 7 is large, then d e C^{flj). We have 



ui + G{x, Dui{x), D^ui{x)) 

>G{x, \Me-^'^^''^Dd{x\XMe-^'^^''^D'^d{x) - X^Me-^'^^''^Dd{x) Dd{x)) 
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in Qj, where p ? is the matrix with entries PiQj. If x G then Xd{x) < 1 and 
so 

(4.12) e"^ < e-^'^'^''^ < 1 in OA- 
Choose M > so that 

(4.13) M{l-e-^) + G{x,0,0)>l in O. 

Then we assume that there is a large A > so that 

G{x, XcDd{x), XcD^d{x) - \^cDd{x) ® Dd{x)) > 
^'^'^^^ for x e fix, Me-^ < c < M. 

Then, by assumption, ui is a classical solution of m + G{x, Du, D^u) > in 0^ and 
Ui = holds (continuously) on dQ. 

To complete the construction, choose C > so that 

(4.15) C + G(a;, 0, 0) > in O 
and 

(4.16) 0<C<M(l-e-i); 

this is possible in view of (4.13). We claim that if u is defined by 



u = 



ui AC in 

C mfl\Qx, 



then it is a supersolution as desired. Indeed, u = on dfl is evident. The constant 
C is a supersolution of the equation in fl by (4.15) and the subset of fl\ on which 
Ui = C does not meet 9(0a) (because of (4.16) and mi = on dQ and Ui = 
M{1 — e~^) when Xd{x) = 1). Since the property of being a supersolution is local 
and closed under finite minimums, we are done. 

We have left the sufficient condition for our construction, (4.14), in a "raw form." 
It is not useful to refine it further without specializing G further. For example, if 
G has the form G{x,p,X) = Go{p,X) — f{x) where / e C{il), in which the x- 
dependence is "separated," the condition that 

lim Go{XcDd{x),XcD'^d{x) - X^cDd{x) ® Dd{x)) = oo 

A— >oo 

uniformly in c in each interval e<c< 1/e, e>0, and x in some neighborhood of 
dfl in fl will certainly allow us to satisfy (4.14). In particular, if Go depends only 
on p, 

lim Gq{p) = 00 

|p|->oo 

and Go(0) — /(a;) < 0, we can uniquely solve the Dirichlet problem for u-\-Go{Du) — 
f{x) = 0. Here we rely on the fact that |£>rf(a;)| = 1 at points of diff'erentiability of 
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d; near d^l (which is all wc need) this follows from the obvious fact that Dd(x) = 
—n{x) on dil where n{x) is the exterior normal to ^1 at a; G dQ. 
In other cases, second-order terms dominate. If 

G{x,p,X) = - trace(^(a;)X) + {b{x),p) - f{x) 

where A, b, f are continuous, and 

{A{x)n{x),n{x)) > on dQ, 

then the left-hand side of (4.14) has the form 

cA^ {A{x)Dd{x), Dd{x)) + 0(A) as A ^ oo. 

Since Dd{x) — —n{x) on dVL, the coefficient of A^ is positive near and it is easy 
to achieve (4.14) by taking A large. Note that we are thus able to assert the unique 
existence of a solution to (4.9) in this case (provided comparison holds) even though 
A may be completely degenerate inside Q,. 

The reader can invent an unlimited array of examples. It should be noted that 
one may often produce both a subsolution and a supersolution by constructions 
like the above and then one need not assume u = is a subsolution. One may also 
usually take maximums and minimums of operators for which (4.14) can be verified 
and stay within this class. Thus, for example, if 

G{x,p,X) = max{-|p|^^-^ - | traceX|^-^ traceX - f{x), - g{x)} 

and < e < 29, a > 0, f,g G C(0), and f,g > 0, then (4.9) has a unique solution. 

Notes on §4. The combination of Perron's method and viscosity solutions was 
introduced by H. Ishii [86]. The definition of a viscosity subsolution (respectively, 
supersolution) u in [86] was that u* is a subsolution (respectively, w,* is a supersolu- 
tion) in the current sense. Solutions are then functions that are both a subsolution 
and a supersolution and continuity is not required. (Note that then the character- 
istic function of the rationals is a solution of u' = 0). With this setup, Perron's 
method does not require the comparison assumption and the statements become 
more elegant (see [86]). 

We mention some other approaches to existence, for even if they are in general 
much more complicated and of a more limited scope, they can be useful in some 
delicate situations. For example, one can use formulas from control and differential 
games to write explicit solutions for approximate equations and then use limiting 
arguments; this approach is used in M. G. Crandall and P. L. Lions [55, Parts III 
and V] and D. Tataru [157]. 

Two other approximations that have been used are discretization and elliptic 
regularization (for first-order equations; this is the method of "vanishing viscosity" 
and its relation to the theory accounts for the term "viscosity solutions"). Having 
solved an approximate problem, one then needs to pass to the limit (with some a 
priori estimates — but see §6!). Existence schemes of this sort have been used in P. 
L. Lions [116], M. G. Crandall and P. L. Lions [53], G. Baxles [14, 15], H. Ishii [82, 
85], P. E. Souganidis [154], and I. Capuzzo-Dolcetta and P. L. Lions [44]. 
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5. Comparison and Perron's method: Variations 

5. A. Comparison with more regularity. Suppose that the supersolution u in 
Theorem 3.3 is more regular; for example, u G C'^{^1) where C^(r2) denotes the set 
of functions that are Holder continuous with exponent 7 G (0, 1] on each compact 
subset of fl. Then the assertions of Theorem 3.3 remain valid if (3.14) is weakened 
to 

, . ( F{y. r. a{x -y),Y)- F{x, r, a{x -y),X) < uj{a\x -y\'^ + \x-y\) 
1 whenever x,y &Ti, r e M, X,Y & S{N), and (3.10) holds 

for some ^ > 2 — 7. Indeed, in the notation of the proof of Theorem 3.3, 

u{x) - v{y) -^\x-yf < u{x) - v{y) - - V? 

implies, upon putting x = y = y and recalling that x will remain in a compact subset 
of n as a ^ 00, a\x — < 2{u{x) ~ u{y)) < C\x — and then a\x — — *■ as 
a 00 provided only > 2 — 7, so the rest of the proof is luichanged. 

More generally and precisely if u is uniformly continuous with the modulus p/2, 
i.e., \u{x) — u{y)\ < p{\x — y\)/2, we have a\x — — p{\:i: — y\) < 0. Defining 
h{a) = sup{r : ar^ — p{r) < 0}, we will succeed if uj{a\x — yf + \x — y\) in (5.1) 
is replaced by uj{ak{\x — y\) + \x — y\) where ak{h{a)) — » as a ^ 00. Since the 
best possible modulus for nonconstant functions behaves like p{r) = cr, in which 
case h{a) = c/a and k{r) = o(r) succeeds, the Lipschitz continuous case is limiting 
with respect to these arguments (however, see below for the case). 

If we assume still more regularity of u, e.g., Du G C^{^1) where < 7 < 1 (with 
7=0 meaning continuity and 7=1 meaning Lipschitz continuity of Du), we need to 
adapt our strategy to obtain further results. 

Let K c Cl he a compact neighborhood of the set of maximum points of u — v 
in O. If Du G there exists a sequence Vn in C{Q) fl C^{K) with the properties 

(5.2) u — ■i/'n— >0 asn— >oo uniformly on fl 
and 

r \Du-Di:n\ < Cn-'^ on K, 

(5.3) I \D^^n\ < Cn^-^ on 

{ \D''Mx)-D^ijn{y)\<Cn^-^\x~y\ ioi x,y G K 

as n ^ 00 where C is a constant. The V'n may then be constructed by mollification 
of u near K and then extended. 
Let {x, y) be a maximum point of 

(5.4) {u{x) - M^)) - {v{y) - My)) - f 1^ - 

as before, we may assume that .r, y € K for large a and n. Since Du{x) — D'il)n{x) = 
a{x — y) and (5.3) holds, we have 

(5.5) a\x — y\< Cn~^ as n ^ oo. 
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Modifying the arguments in the proof of Theorem 3.3 in an obvious way, we will 
still be able to establish comparison if we can show that 



{F{y, r, a{x - y) + D^r,{y), Y + D^^Jniv)) 

- F{x, r, a{x - y) + D^n{x), X + D^iljn{x)))+ 



as n, a ^ oo in some appropriate fashion when r remains bounded and (3.10), 
(5.3), (5.5) hold. 

It is not very informative to try and analyze this condition in general. Instead, 
let us note that in the linear case (see Remark 3.6) 

F{x,r,p,X) = - trace(*S(a;)S(a;)X) + {b{x) , p) + r - f (x) 

with / G C{fl), the first-order terms are harmless from the point of view of verifying 
(5.6) if b is continuous. Via Remark 3.6, the second-order terms contribute at most 

3a trace(*(E(x) - E(y))(S(x) - S(y)) + tTacc{~A{y)D^iPn{y) + A{x)D^iPn{x)) 

when estimating (5.6) above. Here we have put A{x) = *'S{x)T,{x) . Assuming 
E s C^(0) (so A G C^{Q)), we invoke (5.3) and estimate the above expression by 
a constant times 

a\x — y\'^^ + n^~''\x — y\^ + n'^~'^\x — y\. 
By (5.5), this expression may, in turn, be estimated above in the form 

+ n"'"''' ( II as n — > 00. 





Putting a = n^, this becomes 

j^^/3-2A(7+/3) + „-A(7+/3) + l-7 + jl^i^-^^-^^ aS n ^ 00. 

Regarding A as fixed so as to see what we require of 7, we note that we may make all 
of the exponents above vanish by the choices 7 = (1 — 2A)/(1 — A), /3 = 2A/ (1 — A). 

(The case A > i was treated above, so here we are concerned with A < i.) If we 
increase 7 from this value, so that 7 > (1 — 2A)/(1 — A), it follows that comparison 
holds. 

Comparison is thus assured if E G C'^ and Du € under the relation 7 > 
(1 — 2A)/(1 — A). The limit cases A = | and A = need to be discussed separately, 
for 7 = suffices if A = | and 7=1 suffices if A = 0. In the event that 7 = 0, 
the first relation in (5.3) should be replaced by \u — ijjn\ < o(l) and (5.5) by 
C(\x — y\ < 0(1) and then the proof runs as before for 7 = but with the constant 
0(1) appearing everywhere. To treat the other limiting case, A = 0, 7 = 1, note 
that if p is a modulus of continuity of E and A, (5.7) should be replaced by 

- y\f + P{\x - y\) + n\x - y\ 

and (5.5) should read a\x — y\ < C/n. We conclude upon letting n, a — > 00 in a 
manner so that ap{C / {na))^ — > 0. 
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It is worthwhile to compare the uniqueness results above with those that may 
be obtained when solutions are more nearly classical in the sense that they possess 
second derivatives in a strong enough sense. To simplify the exposition, we will 
always assume that (3.13) holds. Then, of course, if we assume that one of u or w is 
(or even everywhere twice differentiablc), the comparison result holds without 
other structure conditions on F beyond properness, since we can work directly with 
maxima of u—v. Unfortunately, this regularity is rarely available, even for uniformly 
elliptic fully nonlinear equations. If we require less regularity of u and v — but still 
much more than assumed in Theorem 3.3 — different structure conditions suffice to 
guarantee the comparison result via finer considerations about the pointwise twice 
differentiability of functions. 

We present below a strategy that establishes comparison without further as- 
sumptions on F if both u + c|a;p and — w + c|a;p are convex for some c > 0. It also 
shows that comparison holds if u, v & ^\oc (^) ^^"^ of u or u lies in W^^^{Q.) 
for some p> N and the structure condition 

r F{x, t,p, X) is uniformly continuous in p 

i uniformly for x G X G S{N), and r,p bounded 

holds. 

In both of these cases, the argument runs as follows: we may assume, without 
loss of generality, that all maxima oi u — v lie in and then general optimization 
results imply that there exists a sequence Pn such that u{x) — v{x) — (pn,x) 
has a strict max;imum at Xn and a;„ x where f is a maximum point of u — v 
over O. Then one knows (of. Lemma A. 3) that if r > is small enough, the set of 
maximum points of u{x) — v{x) — (pn,x) — {q,x) in \x — x„| < r as g ranges over 
the ball Bs contains a set of positive measure if ^ < r~^en where 

£„ = inf [{u{Xn) - v{Xn)) " {Pn, Xn - v) - {u{y) - v{y))] . 
\y-Xn\=r 

Furthermore, either u or u is twice differentiable a.e. (cf. Theorem A. 2) and 

tluTcforc we can find maxima z with [.r,! — z\ < r for some q G Bs that is a point 
where, in the first case, u and v are twice differentiable and in the second case v or 
u is twice differentiable. In the first case we have 

F{z, u{z), Dv{z) +Pn + q, D^v{z)) < < F{z, v{z), Dv{z), D^v{z)). 

In the second case, if v is twice differentiable the same inequality holds and otherwise 
we have 

F{z, u{z), Du{z), D\{z)) < < F{z, v{z), Du{z) -p„-q, D'^u{z)). 

The conclusions are easily reached upon sending 6, r to 0+ and then n ^ oo. 

5.B. Estimates from comparison. We make some simple remarks that hold 

as soon as one has comparison by any means; there are many variants of these 
(including parabolic ones). For example, suppose K > and u — v < K on ^^l 
instead of u < v on dCl. Then u — K is also a subsolution since Jq~^{u — K) = Jq^u 
and F{x,u — K,p,X) < F{x,u,p,X) by properness. Thus u — v < K in Q. In 
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particular, the variant of Theorem 3.3 in which u < v on dCl is dropped and the 
conclusion is changed to u — v < supgQ(u — v)~^ holds. 

In a similar spirit, suppose we have solutions u and v of F{x, u, Du, D'^u) < 0, 
and F{x, u, Du, D^u) > in O where F is proper, satisfies (3.14) and (3.13), K > 0, 

and F{x, r,p, X) + K > F{x, r,p, X). Then w = v + max(supgj2(w — v)^ , K/j) is a 
solution of F{x,w, Dw, D^w) > in Since we have comparison for F, we then 
conclude u — w < supgQ(u — w)~^ = or u — v < max(supgQ(u — v)~^,K/'y). In 
particular, if 

u + G{x,u,Du,D'^u) - f{x) <0 a,ndv + G{x,v,Dv,D'^v)-g{x)>0 in O, 
G is proper and satisfies (3.14), and f,g € C(f2), then 



5.C. Comparison with strict inequalities and without coercivity in u. The 
condition (3.13) was used in the proof of Theorem 3.3 in order to have (3.15). 

If we simply assume there is a 5 > such that cither F(x,u, Du, D'^u) < —S 
or F{x,v,Dv,D'^v) > S we are in the same situation and do not need (3.13) to 
hold. Moreover, if only F{x,u, Du, D'^u) < and F{x,v, Dv, D'^v) > but for 
e > we can find ip^ G C^, > such that IV'el < £ and F{x,u + tJj^,D{u + 
tps), D'^{u + tps)) < —Sg, we conclude that {u + tpg) — v < supgQ(u + tpg — v) and 
then u — v< supgQ('U — v) + 2e and we recover the result as £ | 0. This construction 
can be carried out in some cases. 

5.D. Comparison and existence of unbounded solutions on unbounded 

domains. We first illustrate a method to establish comparison of unbounded solu- 
tions in unbounded domains by showing that if u and v grow at most linearly and 
solve 

(5.8) u + F{Du,D^u) - f{x) <0 and V + F{Dv, D^v) - f{x)>0 in 

where / is uniformly continuous on (i.e., / € UC(M^)), then u< v. After this, 
we will prove that if u is a solution of m + F{Du, D^u) = of linear growth, then 
u e UC(IR^) and finally that Perron's method supplies existence. Thus we will 

prove 

Theorem 5.1. If f £ UC(R^), then u + F{Du,D^u) - f{x) = has a unique 
solution u that grows at most linearly and u e UC(]R^). 

Proof of comparison. The proof proceeds in two steps. First we note that / G 
UC(M^) implies that there is a constant K such that 




(5.9) 



sup {f{x)-f{y)-K\x-y\)<oo 



and then we show that 



(5.10) 



sup {u{x) — v{y) — 2K\x — y\) < oo. 
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By the assumed Unear growth, we have 

(5.11) uix) - viy) < L{1 + \x\ + \y\) on x 

for some L > 0. We choose a family of functions on parameterized by 
i? > 1 with the properties 

f (i) Pr > 0, 

(ii) liminf|j.|^oo/3ij(a;)/|a;| > 2L, 

(iii) \DPr{x)\ + \\D^I3r{x)\\ <Chv R>1, xG M^, 
[ (iv) limfl^^ Pr{x) = for a; G M^, 



(5.12) 



where C is some constant. In view of (5.11) and (5.12)(ii), the function 
*(x, y) = u{x) - v{y) - 2K{1 + \x- y\''Y'^ - {Pr{x) + pR{y)) 

attains its maximum at some point {x,y). Now either (5.10) holds or for large R 

wc have ^{x,y) > 0, which implies 

(5.13) 2K\x-y\<u{x)-v{y). 
Noting that 

(p + Di3r{x), Z + D^PRix)) e J^-+u{x), 
{p - D(3R{y), -Z - D'(3R{y)) e J''-v{y), 

where 

p = {2KD,{1 + |^|2)V2)|^^._., z = {2KDl{l + |^|2)V2)|^^._., 
we see that 

uix) + F{p + D(3r{x), Z + D^0r{x)) < fix), 

v{y) + F{p - D(3R{y), -Z - D^(3R{y)) > f{y). 

From this, using (5.13) and observing that p and Z are bounded independently of 
i? > 1, we have 

u{x) - v{y) < fix) - f{y) + F{p - DPR{y), -Z - D^pR{y)) 

(5.14) - F{p + D(3r{x), Z + D\6n{x)) 

<K\x-y\+C< i(u(x) - v{y)) + C 

where C is a constant independent of i? > 1. Therefore, u{x) — v{y) is bounded 
independently of i? > 1. 

Since ^{x,y) < ^{x,y) < u{x) — v{y), we may then send i? — >■ oo to conclude 
that u{x) - v{y) - 2K{1 + \x- y^)^/^ is bounded and thus that (5.10) holds. With 
the information (5.10) in hand, we repeat the above line of argument, this time 
assuming u{x) — v{x) — 26 > Q holds for some x and putting 

(5.15) Hx,y) = u{x) - v{y) - - y|2 - e{\x\'^ + \y\^) 
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where e, a are positive parameters. For small e we have ^{x,x) > 6 and by (5.10) 

<i> win have a maximum {x, y) at which 

(5.16) 

||i-y|2+£(|x|2 + |y|2) <«(i)-t;(y) <2K|i-y|+C< ^|x-y|2+4ifVa + C 
for some C. Moreover, there are X, Y & S{N) such that 

{a{x -y)+ 2ex, X + 2el) € J^'"^u(f ), {a{x -y)- 2ey, Y - 2el) e J^~v{y) 
and 



(5.17) 




As usual, we have 

u{x) - v{y) < fix) - f{y) + F{a{x -y)- 2ey, Y - 2el) 
^ ' ' -F{{a{x-y) + 2sx,X + 2sI). 

Since 6 < ^{x,x) < u{x) — v{y), and since X <Y hy (5.17), 

6<(jf{\x-y\) + F{a{x - y) - 2sy, X - 2el) - F{a{x -y) + 2sx, X + 2el), 

where w/ is the modulus of continuity of /, and we will obtain a contradiction if 
we show that the iterated limit inferior of the right-hand side as e ]. and then 
a — > cx) is nonpositive. This information is deduced from (5.16), first using the 
extreme inequalities, which imply that a\x — and e(|xp + \y\^) are bounded 
independently of a > 1 and < e < 1. Thus ex,ey and a{x — y) remains 
bounded as e J, 0, while |a; — y| — » as a ^ oo uniformly in £ > 0. Therefore, using 
(5.17) and the uniform continuity of F on bounded sets, 

5 < liminf liminf a;/(|i — ^|) = 0, 

a— >oo ej,0 

which is a contradiction. 

Proof that solutions lie in UC(IR''^). We take the easy way here, using that F is 
independent of x, and note that if ?/ € and u solves u + F{Du, D'^u) — f{x) = 0, 
then w{x) ~ u{x + y) solves w + F{Dw, D'^w) — f{x + y) =0, so that \u{x) — u{x + 
y)\ < sup,gR« \f(z + y)-f{z)\ by 5.B; thus u £ UC(M^). 

Proof of existence. We need only to produce a linearly growing subsolution and 
supersolution and invoke Perron's method. Letting u{x) = A+B{1 + |xp)^/^, we 
note that Du and D^u are bounded so long as B is bounded. By (5.9), |/(a;)| < 
C + K\x\ for some C. If we put B — K and A = supgw \F{Du, D'^u)\ + C, u is a 
supersolution and u = —u is a subsolution. 
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Notes on §5. The discussion of §5. A regarding regularity of Du is new in this 
framework, but such results were first obtained by R. Jensen [102] in a more com- 
plicated presentation. The assertions of §5. A regarding twice differentiable solutions 
(classical or W'^'^, etc.) rely on various versions of the maximum principle as it 
evolved through the works of Aleksandrov [2, 3], Y. Bakelman [6], C. Pucci [143], 
J. M. Bony [37], and P. L. Lions [120]. 

Section 5.C recalls a classical strategy and corresponds to remarks used in M. G. 
Crandall and P. L. Lions [51], H. Ishii [87], H. Ishii and P. L. Lions [96], .... More 
sophisticated uses of estimates from comparison arguments occur in the study of 
regularity questions — sec §10 for references. 

Section 5.D is concerned with the growth of solutions at infinity and its influence 
on comparison- uniqueness results. The relevance of the class UC(R^) for general 
results of this sort was progressively understood in a series of papers by the authors 
and we refer to M. G. Crandall and P. L. Lions [55, Part II] for a few examples 
showing how natural this class is and how it interacts with structure conditions 
on the nonlinearity. Of course, if one restricts the nonlinearity further, other as- 
ymptotic behaviors can be allowed. See H. Ishii [83, 89], M. G. Crandall and P. 
L. Lions [54, 56], and M. G. Crandall, R. Newcomb, and Y. Tomita [59] for more 
information in this direction; in particular, general functions F{x,r,p,X) can be 
used in place of F{p,X) — f{x) (although subtleties concerning existence arise). 

6. Limit operations with viscosity solutions 

Suppose wc have a sequence n = 1,2,..., of subsolutions of an equation 
F = on the locally compact set O. It turns out that the following "limit" 

(6.1) V{z) = Hmsup*'u„(2;) = hm sup{u„(a;) : n> j, x G O, and |-2 — a;| < i} 

in which, roughly speaking, the "limsup" operation and the * operations are per- 
formed simultaneously, will also be a solution of -F < in C Indeed, since U (x) < r 
only if there are e,j > such that Un{z) < r — e for n > j and z € O with 
\z-x\ < 1/j, it is clear that {x € O : U{x) < r} is open in O and thus U e USC(e'). 
We have 

Lemma 6.1. Let Un € USC(C') for n = 1,2, . . . and U be given by (6.1), z & O 
and U{z) < oo. If {p,X) e Jq~^U{z), then there exist sequences 

(6.2) rij -»• 00, Xj e O, {Pj,Xj) e j'^~^Unj{xj) 
such that 

(6.3) {xj,Un^{xj),pj,Xj) {z,U{z),p,X). 

In particular, if each Un is a solution of F < and U < oo on O, then U is a 

solution of F < on O. 

Proof. By definition, there are sequences 

rij — »■ 00, Xj e O such that xj — > z and u„^ {xj) — > U{z); 
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it is also clear that if Zj ^ x in O, then limsupu„^ (^j) < U{x). The result now 
follows at once from Proposition 4.3. 

Remark 6.2. If the assumptions of the lemma are altered to assume instead that 
Un € LSC(O) and {p,X) e J^~v{z), then the assertions are changed by replacing 

(6.4) U.{z) = liminfi,w„(2;) = lim inf{u„(x) : n> j, x £ O, and |-z — a;| < i}, 

n — >oo j — >oo 

putting J^'~ in place of J^'"*" in (6.2) and supersolutions in place of subsolutions. 

Remark 6.3. In fact, the above proof shows more suppose 7i,j is a solution of a 
proper equation F„ < that varies with n. Then the conclusion is that [/ is a 
solution of G < where 

(6.5) G{x,r,p,X) = limintF„(a;,r,p, X); 

note that F„ need not be continuous and that if F„ = F is independent of n but 
discontinuous, then G = F* is the lower semicontinuous envelope of F. Analogous 
remarks hold for supersolutions. 

Remark 6.4. The above results are related to uniform convergence as explained 
next. Let u„ be a sequence of functions on O and 

(6.6) L''(a;) = limsup*M„(x), U_{x) = lim'mi^Un{x). 

n — >oo ^ 

Suppose U{x) = U_{x) on O, lot U(x) denote this common value and assume that 
—00 < U{x) < oo on C Then U is continuous (since it is both upper and lower 
semicontinuous) and lim„^oo Wn(a;) = U{x) uniformly on compact sets. Indeed, if 
this were not the case and uniform convergence failed on some compact set K, there 
would be an £ > and sequences nj oo, Xj G K such that Unj {xj) — U{xj) > s or 
{xj) — U{xj) < —e. Assuming Xj — > x and using the continuity of U , we would 
conclude that \U{x) — ?7(a;)| > £, a contradiction. In order to prove that U = U_ 
one notes that < C/ by definition and typically uses comparison results to prove 
the other inequality. The next result provides a simple example; a typical case in 
which the hypotheses of this result are easily verified is mentioned after the short 
proof. 

Theorem 6.5. Let n be a bounded open set in , H € C(R^), and f € C{Il). 
Consider the problem 

(DP)^ u + H{Du) - eAu = f{x) in ft, u = on dfl, 

and assume that (DP)^ has a subsolution u G C{Q) and a supersolution u G C{Cl) 
independent of e G [0, 1] that vanish on dfl. Then (DP)^ has a solution Ue for 
e G [0, 1] and lim^io "^ei^) — uq{x) uniformly for a; G O. 

Proof. We know that (DP)^ has a solution Ue for e G [0, 1] by Perron's method and 

the assumptions. Moreover, U{x) = \imsnp^^Q*u^{x) and U_{x) = liminfej^o*We(a;) 
are a subsolution and a supersolution of (DP)o by the above and then U < Hhy 
comparison while U > U_hj definition. Thus uq = U = U_ and the convergence is 
uniform. 

Harking back to the discussion of Example 4.6, the reader may easily verify the 
hypotheses of Theorem 6.5 when dfl is of class C"^, H satisfies lim|p|^oo H{p) = oo 
and i?(0) - / < 0. To do so, set m = and u{x) = {\d{x)) A C where A is first 
chosen sufficiently large and then C is chosen in a suitable manner. 
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Notes on §6. The fact that viscosity solutions pass to the limit under uniform 
convergence is an immediate consequence of the definition. This consistency and 
stability property can also be seen as the analogue (in the space of continuous 
functions) of Minty's device for monotone operators as developed by L. C. Evans 
[64, 65]. 

The idea of using only "half-relaxed" limits like U in (6.1) arises naturally when 
attempting to pass to limits with maxima and minima and has been used extensively 
in the calculus of variations and homogenization theory. The use of this concept 
in the area of viscosity solutions was introduced by G. Barles and B. Perthame 
[25, 26] and H. Ishii [90]. This passage to the limit and the associated a posteriori 
uniform convergence through comparison of semicontinuous functions has become 
one of the main features of the theory. Viscosity solutions were first used to prove 
theorems analogous to Theorem 3.2 (but more complicated) with applications by L. 
C. Evans and H. Ishii [68]; they were motivated by questions arising in probability 
(see also §10). Proofs in the spirit we have presented are much simpler. 

7. General and generalized boundary conditions 

7. A. Boundciry conditions in the viscosity sense. In this section we consider 
more general boundary value problems of the form 

r (E) F{x,u,Du,D^u) = in n 
^ ' I (BC) B{x,u,Du,D^u) = on dn 

where (E) denotes the equation and (BC) the boundary condition. Here F is proper, 
as always, and B is a. given function on dfl x M x x S{N) that is also to be proper. 
While it is convenient to allow B to depend on D'^u at this stage, as the reader 
will see, when we turn to existence and uniqueness theorems B = B{x,r,p) will 
be of first order. For example, the Dirichlct condition u — f{x) on dfl arises from 
the choice B{x,r,p, X) = B{x,'r,p) = r — f{x) of B while the Neumann condition 
Un = f{x) (here and later, n{x) denotes the outward unit normal to x £ dfl) 
arises if B{x, r,p) = {n{x),p) — f{x). The Neumann condition is generalized by the 
"oblique derivative" problem, in which B has the form B{x, r,p) = {v{x),p) — f{x) 
where h'{x) is a vector field on dV, satisfying (j^(x), n{x)) > on dft. 

We have already given a meaning to statements like (BC) in Definition 2.2, which 
is, however, clearly inappropriate here since it would only involve the behavior of 
u on dil. The simplest viscosity definition of (BC) we might try is 

Definition 7.1. A function u G USC(^l) is a subsolution of (BC) in the strong 

(viscosity) sense if 

(7.1) B{x,u{x),p,X) <0 forxGdCt,{p,X)GJ^'^u{x) 

and V e LSC(ri) is a supersolution of (BC) in the strong (viscosity) sense if 

(7.2) B{x,v{x),p,X) > for a; e dQ, and {p,X) e J^~v{x). 

Finally, u G C(f2) satisfies (BC) in the strong (viscosity) sense if it is both a 
subsolution and a supersolution in the strong sense. 
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Interpreting (E) in the usual way and (BC) in the strong sense, we obtain notions 
of subsolutions (which will lie in USC(n)), supersolutions (which will lie in LSC(f2)) 
and solutions (which lie in C{i})). 

As indicated by the appearance of the modifier "the strong sense," it will be 
necessary to relax this interpretation to obtain full generality. Before explaining 
this, however, we pause to relate "the strong sense" to the "classical sense." 

It is important to recognize early on that if m G C'^{fl) (i.e., u extends to a twice 
continuously differentiable function in a neighborhood of ft) and 

B{x, u{x), Du{x), D'^u{x)) < for a; G dQ, 

then it is not necessarily true that w is a subsolution of (BC). Indeed, in Remark 

2.7 we determined J'q'^u{x) in the situation where O is an N-submanifold of 
with boundary, x € dil, and dO. is twice differentiable at x. In particular, it follows 
from Remark 2.7 that if 5 < and Xx is given by (2.16) with u in place of tp, then 
{Du{x)—\n{x),X\+ nn{x)'S)n{x)) e Jq'^u{x) for A = /i = and for A > 0, G K. 
Thus, assuming that B is lower semicontinuous, u will not be a subsolution unless 

(7.3) B{x, u{x), Du{x) - Xn{x), Xx + iJ.n{x) n{x)) < for A > 0, e K 

(to see that (7.3) holds with A = and /it arbitrary, just let A J, with /x arbitrary 
but fixed). 

Notice, however, that ifB{x, r,p) is first order and n € M^, then B{x, u{x), Du{x)) < 
implies B{x, u{x), Du{x) — An) < for A > provided 

(7.4) A — > B{x, r,p — An) is nonincreasing in A > 0. 

Moreover, by Remark 2.7, if u is differentiable at x and (p,X) e j'^~^u{x), then 
Du{x) — p is a "generalized normal" to O at x. Since analogous remarks evidently 
hold for supersolutions, we have found a case in which classical satisfaction of (BC) 
implies (BC) in the strong sense. More precisely, we have 

proved 

Proposition 7.2. If B & C{dn x M x R^) satisfies (7.4) for every x e dil, r € 
M, p e andn is a generalized normal toO atx, u & C^(0) andB{x, u{x), Du{x)) < 
{respectively, > 0) forx S dfl, thenu is a subsolution {respectively, supers olution) 
of (BC) in the strong sense. 

This result will be needed in the construction of subsolutions and supersolutions 
below. 

As mentioned above, it will be necessary to relax the interpretation of the bound- 
ary condition "in the strong sense;" one way to motivate this is to consider what 
happens when we try to pass to the limit in the following situation, which is an 
analogue of Remark 6.3 "with boundary:" Suppose F„ G C{n x R x x S{N)), 
Bn € C{dn X R X x S{N)) for n = 1, 2, . . . and w„ be a subsolution of F„ = 
in O, Bn = on dO, in the strong sense. Let F„ — > F, Bn B uniformly on 
compact sets as n — > cxd and u = limsup*^^^ u„ be bounded. What problem does 
u solve? Of course, we know from Remark 6.3 that m is a solution of F < in fi. 
However, when we attempt to check the boundary condition B{x,u{x),p,X) < 

2 -t- 

for X G dCl and {p,X) G u{x), the proof fails. We know, from Proposition 4.3, 
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that (passing to a subsequence if necessary), there exists (p„,X„) S Jjf Un{xn) 
such that {xn,Un{xn),Pn,Xn) — >■ {x,u{x),p, X), but we do not know x„ G dQ (and 
the example below shows this cannot be achieved in general). If a;„ G fi, we have 
Fn{xn,Un{xn),Pn,X„) < (rather than B„(a;„, u„(a;„),p„, X„) < 0). However, in 
this way we can conclude that either 

B{x,u{x),p,X) <0 or F{x,u{x),p,X) <0; 

in other words, u is a subsolution for the boundary condition 

F{x, u, Du, D^u) A B{x, u, Du, D^u) = 

in the strong sense (where a Ab = min{a, b} and a V 6 = max{a, b}). The following 
example, in the context of a Neumann problem, shows that this is the best that 
can be expected. 

Example 7.3. The linear problem 

(7.5) -eu" + u' + u = x + l 

on (0, 1) C M with the homogeneous Neumann conditions 

(7.6) m'(0) = m'(1) = 

has a unique classical solution for £ > and it follows from Proposition 7.2 that 
Ue satisfies (7.6) in the strong sense when we chose the function B that defines the 
boundary condition (7.6) to be given by _B(l,r,p) = p and B{0,r,p) = — p; note 
the monotonicity in the direction of the exterior normals to (0, 1). We will see that 
the limit lim^j^o = u exists uniformly on [0, 1] and the limit u does not satisfy 
(7.6) in the strong sense. 

It is elementary to compute 

Ue[x)=x+ (e^+-e^-)'' X- {e^+ - e^-f " ' 

where A± = (l/2e) (l ± (1 + 4e)'^/'^). Noting that A+ ^ oo and A_ ^ -1 as e i 0, 
one sees that 

Ug{x) u{x) = x-\- uniformly on [0, 1] as e J, 0. 

The function u{x) = x + satisfies 

u' + u = a; + 1 in (0, 1) and u'{0) = 

in the classical sense. However, m'(1) = 1 — 1/e > and therefore u is not a 
subsolution of i?(l,?i'(l)) = in the strong sense. Note, however, that (u' + u — 
(a; + 1)) A w' < does hold at x = 1, in agreement with the above discussion. 

These considerations suggest the appropriate definitions regarding solutions of 
(EVP). 
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Definition 7.4. Let ft be an open subset of R^, F € C{Tl x M x x S{N)), 
and B e C{dQ x R x x S{N)) be proper. Then u is a (viscosity) subsolution 
of (BVP) if ueUSC(n) and 



(7.7) 



F{x,u{x),p,X) < for a; G fl,{p,X) G J^'^u{x), 

F{x,u{x),p,X) A B{x,u{x),p,X) < for a; e dQ, {p,X) G J^'^u{x). 



Similarly, u is a supersolution of (BVP) if m G LSC(O) and 

F{x,u{x),p,X) > for x G ri, (p, X) G J§ u{x), 



(7.8) 



F{x,u{x),p,X)V B{x,u{x),p,X) >0 fov x G dn, {p, X) G u{x). 



'n 

Finally, u is a solution of (BVP) if it is both a subsolution and a supersolution 
Another way to express these definitions is to introduce the functions 



(7.9) G-{x,r,p,X) 
and 

(7.10) G+{x,r,p,X): 



F{x,r,p,X) ifa;GO, 
F{x,r,p,X) hB{x,r,p,X) \{ x G 39. 



Fix,r,p,X) if X G ri, 

F{x,r,p,X)\/ B{x,r,p,X) if x G dn, 



and note that "u is a subsolution of (BVP)" means exactly that u is a solution 
of G- < in ri, etc. Observe that G_ is lower semicontinuous and G+ is upper 
semicontinuous. In fact, G_ is precisely the lower semicontinuous envelope of the 
function 

F{x,r,p,X) if a; Gil, 
B{x,r,p,X) if X G dfl, 



G{x,r,p,X) = 
while G+ is the upper semicontinuous envelope of G. 



7.B Existence and uniqueness for the Neumann problem. We are con- 
cerned here with the following choice B = B{x,r,p) of the boundary condition 
B, 

(7.11) B{x,r,p) = {n{x),p)+fix,r) 

where / G C{d9 x M). With appropriate restrictions on f2, F, and /, we prove that 
(BVP) has a unique solution by a program analogous to that developed in §§3 and 
4. Regarding Cl we will always assume that 

(7.12) O is a compact ./V-submanifold with boundary of M^. 
In addition, we impose the uniform exterior sphere condition 



(7.13) 



3r > such that b{x + rn{x), r) n fl = $ for x G dfl 
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where b{x, r) denotes the closed ball of radius r centered at x. For F we will require 
(3.13), restated here for convenience, 



as well as uniform continuity in (p, X) near d^l; more precisely, we require that 
there be a neighborhood V of dil relative to Q such that 



for some ui : [0, oo) — > [0, oo] satisfying a;(0+) = 0. In addition, we strengthen (3.14) 
to the condition 

(716) i'^(y'^'P'^)r'-^(^'^'^''^) - + + 

^ ' ' \ever x,y&n, r € M, p € M^, X, F G S{N), and (3.10) holds 

where we may use the same w as in (7.15). Observe that (7.16) implies (3.14). For 
/ we assume that 

(7.17) f{x,r) is nondecreasing in r for x £ dO.. 
We have: 

Theorem 7.5. Let (7.11)-(7.17) hold. If u is a suhsolution of (BVP) and v is a 
supersolution of (BVP), then u < v in fl. Moreover, (BVP) has a unique solution. 

The proof consists of verifying the comparison assertion, then producing a sub- 
solution and a supersolution, and invoking Perron's method. We begin with the 

comparison proof. 

Proof of comparison. The strategy involves two types of approximations. As a first 
step, we produce approximations We,Ue such that ^ u, ^ v uniformly as 
s I and Ue,Ve are a subsolution and a supersolution of (BVP) with B{x,r,p) 
replaced by B{x, r,p) + e and B{x, r,p) — e respectively. (Please note that B itself 
will retain its meaning below.) If we then can prove Ue < Ve we find w < t; in the 
limit e I 0. After this, we arrange to deal only with the inequalities involving F 
(and not B) via another perturbation built into the "test functions" we will use. 

Let us proceed with the first step. We need the following elementary lemma 

Lemma 7.6. Let (7.12) hold and v G C(5f},M^) satisfy {n{x),v{x)) > on dn. 
Then there is a function <p e C^(Q) such that 

(7.18) {iy{x),Dip{x)) >1 for x G dn and ip>0 onli. 

The proof of the lemma is deferred until the theorem is established. Choosing 
v = n in Lemma 7.6, let be provided by the lemma; we may assume that the 
support of (p lies in V. Let u be a solution of < in O (that is, a subsolution 
of (BVP)) and v a solution of > 0, and put 

(7.19) Us{x) = u{x) — e(fi{x) — C and Vs{x) = v{x) + e(fi{x) + C 



(7.14) 



^{r-s)<F{x,r,p,X)-F{x,s, p, X) 

for r>s, {x,p, X)eTixR^ X S{N) 



(7.15) 



\F{x, r,p, X) - F{x, r, q,Y)\ < cj{\p -~ q\ + \\X - Y\\) 
for xGV, p,qGR^, X,Y € S{N) 
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where C = C{e) will be chosen later on. Using (7.14) and (7.15), fov x G V and 

2 + 

{p,X) e J[7 Us{x) we have 

F{x,Us{x),p,X) < F{x,u{x),p + eDip{x),X + eD^ip{x)) -^C + w{eM), 

where M = max.-^{\Dip{x)\ + ||D2(/)(a;)||) and, if x G dfl, 

B{x, Ue{x),p) = B{x, u^{x),p + eD(p{x)) — e{Dip(x),n{x)) 
< B{x, u{x),p + eDip{x)) — e. 

Using next that J^y Ue{x) + {eDip{x),£D'^ip{x)) = J-^ u(x) (see Remarks 2.7), we 
conclude that if C = Lo{eM)/^, then Ue is a subsolution of (EVP) with B + £ in 
place of B. Similarly, Ve is a supersolution of (EVP) with B replaced hy B — e. 

We are reduced to proving that if u is a subsolution and is a supersolution of 
(EVP) with B replaced hy B + e and B — e respectively, then u <v. Assume, to 
the contrary, that maxj2(M — w) > 0. We know that maxan(u — v) > maxj2(M — v) 
by §5.B, so there must exist z € dVt such that u{z) — v{z) = maxj2(u — ti) > 0. Put 

^{x,y) = u{x) -v{y) - ^\x - + f{z,u{z)){n{z),x - y) - \x - z\^ on Q x fi, 

and let {xa, Va) be a maximum point of Ey Proposition 3.7, since u(x) — v{x) — 
|x — has a; = ^; as a unique maximum point, 

(7.20) .T„ ^ z, ajxa, — j/ctP ^ 0, u{xa) ^ u{z)^ v{xa) ^ v{z) 
as a ^ cx). For simplicity, we now write (x,y) for (xa,ya) and put 

V'(a;,2/) = -j^\x-y\^ - f{z,u{z)){n{z),x-y) + \x- zf. 

We have not yet invoked the exterior sphere condition (7.13) and do so now. Re- 
stating (7.13) as 

\y — X — rn(x)|^ > for x G 90 and y G O, 

we see it is equivalent to 

(7.21) ('^(a;), y — x) < ■^\y — for x G dO. and y G il. 

Using (7.21) we compute that if f G dO. then 

B{x,u{x),Dxip{x,y)) = B{x,u{x),a{x - y) 

-/(z,u(z))n(z)+4|x-zp(x-z)) 

> --^\x - - f{z, u{z)){n{x), n{z)) 

+ 0(|x-zn+/(x,u(x)) 

and in view of (7.20) this implies 

B{x, u{x), Dxtp{x, y)) > o(l) as a ^ oo if x G dCl. 
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A similar computation (this time involving (7.17)) shows that 

B{y, v{y), -Dyipix, y)) < o(l) as a -»• oo if y G dCl. 

Recalling that u is a subsolution with respect to B + e, etc., and the definition of 
subsolutions and supersolutions, we conclude 



(7.22) 



F{x, u{x), D^tl){x, y),X)<0 for {D^'ip{x, y),X) e J§"^u(f ), 



F{y, viy), -Dy^ix, y),Y)>0 for (-D^^(f , y),Y) & J§ v{y), 



provided a is large. 

The remaining arguments are very close to those given in the proof of Theorem 
3.3. We apply Theorem 3.2 with k = 2, Oi = O2 = ^, Ui = u, U2 = ~v, and 
£ = 1/a to find X,Y G S{N) such that 



{D^ipix, y),X)e J§+, {-Dy^{x, y),Y)& v{y) 



I 0\ ^ ( X \ ^ , , 1 ,2 



and 



(7.23) -(a+p||)(^^ l)^[o -VJ^^+a^ 

where A = D'^tp{x,y). Observing that 

^ = "(-7 7^) = 2c^' (i/ ~/ )+0{a\x-z\' + \x-z\''), 
we see from (7.23) that 

-3a ^^+o(l)<^^ _^Y^-^^(^-j as a ^ 00. 

We want to use this information in conjunction with (7.22); to this end, we formulate 

the next lemma (which will be used later as well). 

Lemma 7.7. Let (7.15) and (7.16) hold. Let 6 > 0, x,y e V , r & R, p & , 
X,Y gS{N), and 



,0 /;-vo -Y J - \ -I I J ' \o I 

hold. Then 

F{y, r,p, Y) - F{x, r,p, X) < u;{{a + |5)|a; - y\^ + \x- y\i\p\ + 1)) + 2uj{5). 
Proof. By assumption, 
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Using (7.16) we therefore have 

F{y,r,p, Y + SI)- F{x,r,p, X - SI) < a;((a + p)\x - + \x- y\{\p\ + 1)). 

Next, (7.15) implies 

F{y,r,p,Y + SI) - F{x,r,p,X - 61) > F{y,r,p,Y) - F{x,r,p, X) - 2uj{S) 

and the result follows. 

We complete the proof of comparison. Prom the lemma and the developments 
above, 

< F{y, v{y), -Dy^{x, y),Y) - F{x, u{x), D,^{x, y),X) 

< F{y, u{x), -DyiPix, y),Y)- F(x, u{x), -Dy^{x, y), X) 
- -f{u{x) - v{yj) + w(4|x - zp) 

< u;{a\x - + \x- y\{\Dy^l,{x, y)\ + 1)) - j{u{x) - v{y)) + o(l) 

as a oo. From this wc obtain a contradiction as in the proof of Theorem 3.3. 

Proof of existence. In order to establish existence, we may use Perron's Method 
from Theorem 4.1 as supplemented by Remark 4.5 to deduce that it suffices to 
produce a subsolution and a supersolution. Using Proposition 7.2, we sec it suffices 
to produce a classical subsolution and supersolution. Let (fi e C^(f2) be as in 
Lemma 7.6. Put 

u{x) = -Aip{x) - B, u{x) = Alp{x) + B 
where A, B are constants to be determined. Observing that 

B{x,u{x),Du{x)) > A{n{x),Dip{x)) + f{x,0) >A + f{x,0) 
for X G dfl and 

F{x,u{x),Du{x),D^u{x)) > F{x,0, ADip{x), AD'^ifix)) + 'yB 
for X &fl we see that if we first choose A and then B by 

A = sup|/(x,0)|, B=-_ sup \F{x,0,p,X\, 

9^ 'y n,|p| + ||X||<AM 

where M = su\yQ{\D(p{x) + \\D'^(p{x)\\), then u is a supersolution. Likewise, u will 
be a subsolution. It remains to prove Lemma 7.6. 

Proof of Lemma 7.6. In view of (7.12), there is a function ip on such that 
> on M''^\r2, < on n, and \Dip\ > on Multiplying by a constant 
if necessary, we may also assume that \Dtp{x)\ > {u{x),n{x))~^ for x G dQ. By 
standard approximation theorems, we see that there is a sequence of functions 
ipk such that i/'fe ^ "0 in C^(fl) as fc ^ oo. In particular, we have 

{i'{x),Dipk{x)) {u{x),Dtjj{x)) uniformly on dil as A; — > oo. 

Since n{x) = D'tjj{x) /\Dijj{x)\ and \D'ilj{x)\{i'{x),n{x)) > 1 for x e dfl, if k and 
then C are sufficiently large, (fi = tpk + C has the desired properties. 
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7.C. The generedized Dirichlet problem. We turn to the Dirichlet problem 

( F(x,u,Du,D'^u) = inO, 
{ u — j{x) =0 on as 2 

where / e C{dn). 

The simplest example of this problem with a unique solution that does not 
satisfy boundary conditions in the strong sense is the following: F{x, u, Du, D^u) = 
u — h{x) and f{x) =0. If /i is continuous on fi, then u = h ]s evidently the only 
solution of this problem and it does not satisfy u = on dO. in the strong sense 
unless /i = on dil. 

We will provide a comparison theorem, Theorem 7.9 below, with this example 

as a special case. This comparison theorem, however will be different in character 
from those we have obtained so far in that it will assert comparison for a continuous 
subsolution and supersolution. As we know, the ability to compare semicontinuous 
subsolutions and supersolutions allows us to prove existence. Thus the following 
counterexample to existence indicates that the restriction to continuous functions 
(or some other assumption) is necessary in Theorem 7.9. 

Example 7.8. Let iV = 2 and Vl = {{x,y) : -1 < a; < 1, < y < 1}. Let 
/ e C{dn) satisfy < / < 1, /(a;,0) = 0, and /(x, 1) = 1 for \x\ < 1. We claim 
that the problem 



, + xUy = in n, 

7-25 / «o 

J on oil 



{" 



does not have a solution. To see this, suppose the contrary and let u{x^ y) be a 
solution. We claim that then for fixed |a;| < 1, v{y) = u{x, y) is a solution of 

(7.26) v + xv' = ()m{(),l), w(0) = 0, v{l)^l. 

Assume for the moment that this is the case. Then we remark that Theorem 7.9 
below shows that v is uniquely determined by (7.26), i; = is a solution of (7.26) 
if x > and v{y) = e^^~y'>/^' is a solution of (7.26) if x < (these last assertions 
the reader can check by computations). Thus we see that u(a;, y) = if x > and 
u{x,y) = e^-^"")/^ if a; < 0; however, u cannot then be continuous at (0, 1) and we 
conclude that (7.25) does not have a solution. 

It remains to verify that if u is a solution of (7.25), then w is a solution of 
(7.26). One way to check this is as follows: Fix x & (—1)1) and set v{y) = 
u{x,y). Let (p G C'^dO, 1]) and assume that v — ip has a unique maximum at 
some point y S [0, 1]. For a > let {x,y) be a maximum point of the function 
u{x, y) — (f{y) — a\x — xp. It follows from Proposition 3.7 that then x ^x, y -^y 
as a ^ oo; moreover, we have u{x,y) + x(p'{y) < 0. Prom this information we 
deduce that v{y) + xp'{y) < and conclude that v is a subsolution of (7.26) with 
X = X. In a similar way, one shows that v is a supersolution. Finally we remark 
that u = 1, M = are, respectively, a supersolution and a subsolution of (7.25). 
Thus comparison of semicontinuous scmisolutions of (7.25) must fail. (If the reader 
is concerned about regularity questions and the "corners," note that the corners 
can be smoothed up without changing the essential points of the discussion.) 

We turn to the comparison result for continuous solutions. 
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Theorem 7.9. Let (7.12), (7.14), (7.15), and (7.16) hold. If u, v G C{n), u is a 
subsolution of (7.24) and v is a supersolution of (7.24), then u < v. 

Proof. As before, we argue by contradiction and so suppose that ma,x^{u — v) >0. 
We may assume that ma.x-^{u — v) = u{z) — v{z) > for some z G dCl. We divide 
our considerations into two cases. 

First, we consider the case when v{z) < f{z). For a > 1 and < e < 1, we 
define the function 

$(a;, y) = u{x) — v{y) — \a{x — y) + en{z)'\^ — e\y — z\^ on f2 x Q. 

Let (x, y) be a maximum point of $. We assume a is so large that z—(e/a)n{z) € Q.. 
The inequality $(i,y) > $(z — {e/a)n{z),z) reads 

\a{x — y)+ sn{z)f' + e\y — zf' < u{x) — v{y) —u(^z— —n{zfj + v{z) 

from which it follows (using the continuity of u\) that if e is fixed, x,y ^ z and 
a{x — y)+ en{z) — >■ as a ^ oo. Indeed, it is clear that a{x — y) remains bounded 
as a — *■ oo, so assuming that a{x — y) ^ w, x,y ^ z (along a subsequence) as 
a — > 00, one finds 

\w + sn{z)f + e\z - zf < u{z) - v{z) - u{z) + v{z) < 0, 

where the last inequality follows from the definition of z, whence the claim is im- 
mediate. From this we see that x = y~ {£n{z) + o(l))/a as a — > oo and hence that 
X € Cl if a is large enough. Since v{z) < f{z), we have (using the continuity of v) 
v{y) < f{y) if a is large enough. Thus, if a is large enough, we have 

F{x,u{x),p,X) <0 for (p,X) €J^"^u(x) 

and 

F{y, v{y), q,Y)>0 for (g, Y) G J^^~v{y). 
Observe that if we set (p{x, y) = \a{x — y) + en{z)\'^ + e\y — z\'^, then 

D^(p{x, y) = 2a{a{x - y) + en{z)), 

-Dyip{x, y) = 2a{a{x - y) + £n{z)) - 2e{y - z), 

and 

Z)V(x,y)=2a^(_^^ 'l)+2e(l . 

Now, using Theorem 3.2 together with Lemma 7.7, calculating as usual, sending 
a — > 00 and then £ J, 0, we obtain a contradiction. 

It remains to treat the case when v{z) > f{z). Since u{z) > v{z) this entails 
u{z) > f{z). Replacing the above $ by 

<^{x,y) = u{x) — v{y) — \a{x — y) — en{z)\'^ — e\x — z]"^, 

we argue as above and obtain a contradiction. 
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7.C'. The state constraints problem. The problem 



(7.27) 



F{x,u,Du,D^u) < in O, 
F{x, u, Du, D'^u) > in n 



corresponds to an important problem in optimal control, the so-called state con- 
straints problem provided F has the form of the left member of (1.10). This is, 
indeed, the extreme form of (7.24) where f{x) = — oo. To make this clear, we give 
the definitions precisely. A function u G USC(ri) (respectively, v G LSC(n)) is 
called a subsolution (supersolution) of (7.27) if 



2 -t- 

F{x,u(x),p,X) < for a; G and (p,X) G J ' u{x) 



(respectively, 

F{x,v{x),p,X) > for a; G and {p,X) €l^~v{x)). 
An application of Theorem 7.9 immediately yields 

Theorem 7.10. Let (7.12), (7.14), (7.15), and (7.16) hold. Ifu,ve C(II), u is a 
subsolution of (7.27) and v is a supersolution of (7.27), then u < v on Q. 

Indeed, u is a supersolution and w is a subsolution of (7.24) with f{x) > u{x) on 

In fact, a subsolution v of (7.27) is a supersolution of (7.24) for any / £ C{dil) 
and hence a "universal boimd" on all solutions. This implies that solutions of (7.27) 
do not exist in general. For example, if _F is a linear uniformly elliptic operator, 
(7.24) is uniquely solvable for arbitrary / G C{dil) and then there can be no uni- 
versal bound that is continuous in il. On the other hand, if F is suitably degenerate 
in the normal direction at dfl, one can sometimes find such a supersolution. 

7.D. A remark (BC) in the classical sense. The reader will have noticed by 
now that we gave some examples of solutions of (BC) that were smooth but did 

not satisfy the boundary condition in the classical sense. In these examples, the 
equation (E) is rather degenerate. In fact, "degeneracy" of a sort is necessary for 
this to happen, as we now establish. 

Proposition 7.11. Let B G C{dQ. x M x R^), F he proper, and u G £"^($1) he a 
subsolution [respectively, supersolution) of (BVP) . Lf 



(7.28) 

[respectively, 
(7.29) 



lim sup^^oQ F[x,r,p,X — jin ^ n) > 

for [x,r,p,X) GdnxRxR^ X S[N),n € R^\{0} 



lim inf ^^oo F[x, r,p,X + fin (g) n) < 

for [x,r,p,X) ednxRxR^ X S[N),nG M^\{0}), 



then u is a classical subsolution [respectively, supersolution) of [BC). 
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Proof. It will suffice to treat the subsolution case. We rely again on Remarks 
2.7(iii). Let y e M^\f2 and z e 5f2 be a nearest point to y in Q. It is clear 
that the set of such nearest points (as y varies) is dense in dQ. Then Cl c {x G 
: — y| > — ^;|} = O where the last equality is the definition of O; now 

this inclusion evidently implies Jq Zero C J-^ Zero(z). On the other hand, we 
computed J^'^Zero(z) in Remark 2.7 — with the notations n = [y — z)/\y — z\ and 
r = \y — z\, this computation shows that 

^— An, — / + /xn <8) S jQ'''Zero(2:) 

for A > and G M. Thus the assumption that « is a supersolution implies that 

B{z, u{z), Du{z) - An) A F ^z, u{z), Du{z) - An, D'^u{z) + ^/ + /in n^ < 

for A > and /x G R. Taking the limit superior as ^ — c» and using the fact that 
assumption (7.28) implies B{z,u{z), Du{z) — An) < for A > and then letting 
A I wc find B{z,u{z), Du{z)) < 0. Since z was an arbitrary "nearest point" and 
these are dense in dQ, the proof is complete. The reader will notice that it suffices 
to choose n in (7.28) and (7.29) from N-^{x). 

7.E. Evilly nonlinear boundary conditions. We conclude this section by de- 
scribing an extremely general existence and uniqueness result for a fully nonlinear 
ffist-order boundary operator of the form B = B{x,p), B S C(dQ, x M^). We will 

require that B satisfies 

(7.30) B is uniformly continuous in p uniformly in x G dO, 
and 

(7.31) \B{x,p) - B{y,p)\ < uj{\x - y\{l + \p\)) 

for some lo : [0, oo) —>■ [0, oo) satisfying lu{0+) — and for some > 

(7.32) B{x,p + Xn{x)) > B{x,p) + vX for x&dVl, A > 0,p e K^. 
We will also need to strengthen a bit the regularity of O by assuming that 

(7.33) 90 is of class C^'^ 

Observe that (7.33) is stronger than (7.12) and (7.13) since is bounded. The 
loss of generality, however, is not really more restrictive for applications, since O 

typically is cither smooth or has "corners," and corners arc not allowed by (7.12) 
in any case. (Corners typically require some ad hoc analysis.) 
We have 

Theorem 7.12. Let (7.30)-(7.33) and (7.14)-(7.16) hold. If u is a subsolution of 
(EVP) and v is a supersolution of (EVP), then u < v in Q. Moreover, (EVP) has 
a unique solution. 
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Example 7.13. The Neumann boundary conditions of (7.11) arc special cases of 
those considered above; note that (7.31) holds since n is Lipschitz continuous on dO, 
because of (7.33). Similarly, this setting includes the more general linear oblique 
derivative boundary condition corresponding to 

(7.34) B{x,p) ^ {j{x),p) + fix) 
where 7 is Lipschitz and satisfies 

(7.35) {-f{x),n{x)) >u>0 for x G dCt; 

note that / is merely required to be continuous on dQ. 

Nonlinear examples arising in optimal control and differential games are typically 
written in the form 

(7.36) B{x,p) = sup{{j„{x),p) + U{x)) 
and 

(7.37) B{x,p) = inf sup((7„,/3(a;),p) + faA^)) 

P a 

where 7^, fa, 7a;,/3, fa.fj sue functions that satisfy the conditions (including conti- 
nuity) laid on 7, / above uniformly in the parameters a, f3. 

As a last example, we offer the famous capillarity boundary condition that cor- 
responds to 

(7.38) B{x,p) = {n{x),p) - k{x){1 + \pf)'/' 

where k S C{dfl) and sup^Q I'^l < 1- 

We will not prove Theorem 7.12 since its proof, while based on the strategies 
presented above, requires some rather complex adaptations of them. We remark 
that the form of the boundary operator B{x,p) could be generalized to B{x,r,p) 
above, provided that B is nondecreasing with respect to r G R and (7.30)-(7.32) 
hold uniformly in bounded r. 

Notes on §7. Boundary conditions in the viscosity sense first appeared in P. 
L. Lions [121] and B. Pertliamc and R. Sanders [142] for Nemnann or oblique 
derivative problems for first order equations, making this problem well posed! A 
general formulation was given later by H. Ishii [90]. Existence and uniqueness 
was shown in P. L. Lions [121] (for the first-order case), H. Ishii and P. L. Lions 
[96] (second-order case), and P. Dupuis and H. Ishii [60, 61] (second-order and 
nonsmooth domains). However, it is worth emphasizing that the presentation we 
give using Theorem 3.2 simplifies the arguments substantially (even in this special 
case). 

Concerning the Dirichlet problem, most of the works required (and often still 
require) continuous solutions up to the boundary and prescribed data on the entire 
boundary. It was first noted in P. L. Lions [116], G. Barles [14, 15] and H. Ishii 
[86] that this can be achieved for special classes of equations by imposing compat- 
ibility conditions on the boundary data or, in general, by assuming the existence 
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of appropriate super and subsolutions. It was also noted in M. G. Crandall and R. 
Newcomb [58] and P. E. Souganidis [156] that some part of the boundary may be 
irrelevant (in the case of a first-order equation). 

Progress on the understanding of Dirichlet boundary conditions was stimulated 
by the "state constraints" problems studied first by M. Soner [149, 150] and later 
by I. Capuzzo-Dolcetta and P. L. Lions [44]. This led to the "true" viscosity for- 
mulation of Dirichlet conditions as considered in H. Ishii [90] and G. Barles and 
B. Perthame [25, 26]. This formulation, in some sense, automatically selects the 
relevant part of the boundary for degenerate problems and yields uniqueness for 
continuous solutions. However, the existence of a continuous solution fails in gen- 
eral (Example 7.8), and the situation is not entirely clear except for first-order 
optimal control problems [25, 26, 90]. 

We also mention at this stage that some uniqueness results for semicontinuous 
solutions have begun to emerge (G. Barles and B. Perthame [25, 26] and E. N. 
Barron and R. Jensen [33, 34]) for optimal control problems and the study of some 
second-order nondegenerate state constraints problems (J. M. Lasry and P. L. Lions 
[111]). 

Finally, the possibility of solving general equations with general fully nonlinear 
oblique derivative type boundary conditions a rather startling fact was illus- 
trated in G. Barles and P. L. Lions [24] for first-order equations. The full second- 
order result of Theorem 7.12 is a generalization by G. Barles [21] of a similar result 
of H. Ishii [91]. 

8. Parabolic problems 

In this section we indicate how to extend the results of the preceeding sections 
to problems involving the parabolic equation 

(PE) ut + F{t,x,u,Du,D'^u)=0 

where now u is to be a function of {t, x) and Du, D^u mean Dxu{t, x) and Dlu{t, x). 
We do this by discussing comparison for the Cauchy-Dirichlet problem on a bounded 
domain; it will then be clear how to modify other proofs as well. Let O be a locally 
compact subset of K^, T > 0, and Ot = (0,T) x O. We denote by + , P'^~ 
the "parabolic" variants of the semijets Jq~^, Jq ', for example, if u: Ot —>■ R then 
P^'+u is defined by {a,p,X) e M x x S{N) lies in ■p^^u{s,z) if {s,z) G Ot 
and 

u{t, x) <u{s, z) + a{t — s) + (p, X — z) + ^{X{x — z),x — z) 
^^'^^ + o{\t - s\ + \x - z\^) a.s Ot 3 {t,x) ^ is,z); 

similarly, Vq~u = —'Pq'^{—u). The corresponding definitions of , V^q are 
then clear. 

A suhsolution of (PE) on Ot is a function u G USC(Ot) such that 
(8.2) a + F{t,x,u{t,x),p,X) < for {t,x) e Ot and ia,p,X) e 'P'^^u{t,x); 
likewise, a supersolution is a function v € LSC(Ot) such that 



(8.3) a + F{t,x,v{t,x),p,X) >0 for {t,x) e Ot and {a,p,X) eV^ v{t,x); 
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and a solution is a function that is simultaneously a subsolution and a supersolution. 

Remark 8.1. Suppose we set xo = t and x = {xq, x). Then (PE) is an equation of the 

form F{x, u, D^u, D'iu) = u^o "I that is proper if and only if F{t, x, u, Du, D'^u) 

is proper when t is held fixed. The definition of subsolutions, etc., of (PE), which 
takes into account that only the first derivative with respect to t appears, does not 
coincide with the definition of subsolutions, etc., of = under this correspon- 
dence. However, it is not hard to see that the two notions are equivalent. 

We will illustrate the additional considerations that arise in dealing with the 
Cauchy-Dirichlet problem for (PE) as opposed to the pure Dirichlet problem in the 
elliptic case. The problem of interest has the form 

' (E) ut + F{t, X, u, Du, D\) = in (0, T) x fl, 

(8.4) I (BC) u{t, a;) = for < i < T and a; e dfl, 

, (IC) u{0, x) = tp{x) for a; G n, 

where C is open and T > and -0 G C(0) arc given. By a subsolution 
of (8.4) on [O.T) x Q we mean a function u G USC([0, T) x fi) such that u is a 
subsolution of (E), u{t,x) < for < f < T and x G dil and u{0,x) < i^{x) for 
X G O-the appropriate notions of supersolutions and solutions are then obvious. 

Theorem 8.2. Let Vl dM.^ he open and hounded. Let F G C([0,T] x 17 x R x 
x S{N)) he continuous, proper, and satisfy (3.14) for each fixed t G [0,T), with 
the same function . Ifu is a suhsolution of {8 A) andv is a supersolution o/(8.4), 
then u <v on [0,T) x CI. 

To continue, we require the parabolic analogue of Theorem 3.2. It takes the 

following form 

Theorem 8.3. Let Ui G USC((0,T) x Oi) for i = 1, . . . ,k where Oi is a locally 
compact suhset of M^* . Let he defined on an open neighborhood of (0, T) x 

Oi X • • • X O}. and such that {t,xi, - ■ ■ , x^.) (pit, Xi, ■ ■ ■ , x^) is once continuously 
dijferentiahle in t and twice continuously differentiahle in (xi, • • • , x^) G Oi x ■ • ■ x 
Ok- Suppose that i G (0, T), f , G Oi for i = 1,. . . ,k and 

w{t,xi, ...,Xk)= ui{t,xi) H \-Uk{t,Xk) - ip{t,xi, ...,Xk) 

< w{i,xi, ...,Xk) 

for < t <T and Xi € O. Assume, moreover, that there is an r > such that for 
every M > there is a C such that for i = 1, . . . ,k 

hi <C whenever (bi,qi,Xi) e P'^^u-i_{t. Xi). 

\xi - Xi\ + \t-i\<r and \u,{t,x,)\ + \q,\ + \\Xi\\ < M. 

Then for each e > there are Xi G S{Ni) such that 

2 + ^ 

(i) {bi,Dxiip(t,xi,. . . ,Xk),Xi) G Va^ u^{t,Xi) fori = l,... ,k, 



(8.6) 



(ii) -(^ + ||A||V< f ; ••• M <A + eA' 



e 



... Xk. 

^ (iii) 6i H \-bk= V>t{i,xi, . . .,Xk), 

where A = {Dlip){i,xi, . . .,Xk). 
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Observe that the condition (8.5) is guaranteed by having each Ui be a subsolution 

of a paraboHc equation. 

Proof of Theorem 8.2. We first observe that foY e > 0, u — u — e / {T ~ t) is also a 
subsolution of (8.4) and satisfies (PE) with a strict inequality; in fact, 

Du,D^u)< 



Since u <v follows from u < v in the limit £ J, 0, it will simply suffice to prove the 
comparison under the additional assumptions 



(8.7) 



(i) Ut + F{t, X, u, Du, D'^u) < -e/T^ < and 

(ii) \miu{t,x) = — oo uniformly on CI. 



We will assume 

(8.8) {s,z)e{0,T)xn and u{s, z) - v{s, z) ^ S > 

and then contradict this assumption. We may assume that u, —v are bounded 
above. Let {t,x,y) be a maximum point of u{t,x) — v{t,y) — {a/2)\x — j/p over 
[0, T) X O X O where a > 0; such a maximum exists in view of the assumed bound 
above on u, —v, the compactness of fl, and (8.7)(ii). The purpose of the term 
{a/2)\x — 2/p is as in the elliptic case. Set 

(8.9) M„ = u{i, x) - v{i, y)-^\x- y\\ 
By (8.8), Mc > 5. If i = 0, we have 

< (5 < M„ < sup U{x) - V'(y) - f k - y\A ; 

however, the right-hand side above tends to zero as a — > c» by Lemma 3.1, so f > 

if a is large. Likewise, x,y E D, ii a is large by u < v on [0, T) x dfl. Thus we may 
apply Theorem 8.3 at {i, x, y) to learn that there are numbers a, b and X,Y E S{N) 
such that 

2 + ^ 2 — 

{a,a{x - y),X) eVq u{t,x), {b, a{x - y),Y) e Vq v{t,y) 

such that 

(8.10) a -6 = and -^ct(j^ < (^^ ~/ 

The relations 

a + F{i, X, u{t, x), a{x — y), X) < — c, 
b + F{t,y,v{i,y),a{x - y),Y) > 0, 

and (8.10) imply 

c < F{i, y, v{i, y), a{x - y), Y) - F{t, x, u{t, x), a{x -y),X) 
<uj{a\x-y\ + \x-y\) 
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which leads to a contradiction as in the proof of Theorem 3.3. 

Let us mention a couple of other adaptations of results above to parabolic prob- 
lems. Section 5.B may be regarded as establishing continuity of solutions with 
respect to boundary data and the equation itself. In the parabolic case, there is the 
initial data, the boundary data, and the equation to consider. In the above context, 
we may consider a solution of Ut + F{t, x, Du, D'^u) < in (0, T) x f2 and a solution 
of Vt + G{t,x,Dv,D^v) > in (0,T) x with the continuity and boundedness 
properties assumed in Theorem 8.2. Suppose g{t) > {G{t, x,p, X) — F{t, x,p, X))+ 
where g is continuous and {u{t,x) — v{t,x))~^ < Ki for {t,x) S (0,T) x 90, 
(u(0,a;) — w(0,a;))+ < K2 for x & Q.. Then the function w{x,t) = v{x,t) + 
vaax{Ki, K2) + /J g{s) ds is a solution of Wt + F(t, x, Dw, D^w) > 0. In this way, 
if we have comparison, we conclude that 

u{t,x) <v{t,x) +mayi I sup (u — v)~^ , sup (u — v)~^ \ + / g{s)ds. 

\(o,T)xdn (o\xa J Jo 



In many cases we may put 



g{t)= sup {G{t,-)-F{t,-))+, 

and a simple example is G{t, x,p, X) = F{t, x,p, X) + f{t, x). 

Notes on §8. The presentation follows M. G. Crandall and H. Ishii [48]. Let 
us also mention that it is well recognized that most results concerning stationary 
equations have straightforward parabolic analogues; indeed, the parabolic situation 
is often better, since 

Ut + F{x, u, Du, D'^u) = 

more or less corresponds to Ait + F{x, u, Du, D'^u) = with large A > 0. Moreover, 
the special linear dependence on ut in the parabolic case allows one to let F depend 
on t in a merely measurable manner; see, e.g., H. Ishii [84], B. Perthame and P. L. 
Lions [129], N. Barron and R. Jensen [32], and D. Nunziante [137, 138]. Another 
special feature of the parabolic case is that, owing to the special structure, one 
can often allow rather singular initial data. For example, infinite values may be 
allowed and semicontinuity may suffice, depending on the situation (see, e.g., M. 
G. Grandall, P. L. Lions, and P. E. Souganidis [57] and E. N. Barron and R. Jensen 
[33, 34]). 

9. Singular equations: an example from geometry 

Let p® q — {piqj}, the matrix with entries PiQj- It can be shown that if i/i is a 
smooth function and Dijj does not vanish on the level set F = {'0 = c} and u is a 
classical solution of the Cauchy problem 

(9.1) Ut - trace (^(^I - D^u^ = 0, u{0,x) = ipix) 

on some strip (0, T) x R^, then Ft = {u{t, •) = c} represents the result of evolving F 
according to its mean curvature to the time t, whence there is geometrical interest 
in (9.1). Indeed, in less regular situations where we have viscosity solutions, has 
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been proposed as a definition of the result of evolving T in this way. The nonlinearity 
involved, which we hereafter denote F{p,X) = — trace((/ — {p ^ p)/\p\^)X), is 
degenerate elliptic on the set p ^ and undefined at p = 0. Thus the results of 
the preceeding sections do not apply immediately. However, this is easy to remedy 
using the special form of the equation. The extensions of F to (0, X) given by 

. . ( F{p,X) if p 7^0, , ( F{p,X) ifp^O, 

9.2 F(p,X)^l „ ' F(p,X)^l ''\, ' ' 

^ ' ' 1 -2||X|| ifp = 0, ^ ' \2\\X\\ ifp = 

are lower semicontinuous and upper semicontimious respectively. We define u to 
be a subsolution (respectively, supersolution) of ut + F{Du, D^u) = if it is a 
subsolution of Ut + F_{Du, D'^u) = (respectively, of Ut + F{Du, D^u) = 0) and a 
solution if it is both a subsolution and a supersolution. The reason that this will 
succeed is roughly that, in the analysis, we will need only to insert X = when we 
have to deal with p = 0, so what one does with p = is not important so long as 
it is consistent. To illustrate matters in a slightly simpler setting, let us consider 
instead the stationary problem 



(9.3) u + F{Du,D\)-f{x)=0 in 



where we use the corresponding definitions of subsolutions, etc. Theorem 5.1 re- 
mains valid for the current F. 

Theorem 9.1. Let f e UC(M^). Then ((9.3) has a unique solution u G UC(M^). 

Wc sketch the; proof, which proceeds according to the outline given in the proof 

of Theorem 5.1 with slight twists. 

Proof of comparison. The comparison proof is a slight modification of that of §5.D. 
We begin assuming that u,v are a subsolution and a supersolution of (9.3) and 
u{x) — v{y) < L{1 + \x\ + |y|) (which is (5.11), and then proceed as before, ending 
up with 

u{x) - v{y) < fix) - f{y) + F{p - Dpniv), -Z - D^pniv)) 

-F{p + DPr{x),Z + D^(3r{x)) 

in place of (5.14). This still implies a bound on u{x) — v{y) — 2K\x — y\ as before. 
Then proceeding still further with the proof, we adapt slightly and consider a 
maximum point (i, y) of 

^{x, y) = u{x) - v{y) - {a\x - yf + e{\x\^ + \y\^)), 

which will exist by virtue of the bound already obtained. We assume, without loss 
of generality for what follows, that ^{x, y) > so that for some C 

(9.5) a\x - + e(|x|^ + \yf) < u{x) - v{y) < 2K\x -y\ + C. 
Using Remark 3.8, we then have the existence of X, F G S{N) such that 

(9.6) {p + 2ex, X + 2el) e J^'"^u(5), (p - 2ey, Y - 2sl) e J^'~v{y) 



54 



M. G. CRANDALL, HITOSHI ISHII, AND PIERRE-LOUIS LIONS 



where 

P = ^a\x-y\'^{x-y), 
||^||,||1"|| < CiQ!|f and X <Y. 

Now (9.5) implies that s{\x\^ + \y\^) and a\x — y\^, and hence p, are bounded 
independently of £ < 1 for fixed a > 1 while — y]^ is also bounded independently 
of a. Hence, ex,ey as £ J. 0. We have the following analogue of (5.18) 

u{x) - v{y) < (Fip - 2sy, Y - 2el) - F{p + 2ex, X + 2sl)) + f{y) - f{x) 

that leads, by use of the estimates above, to 

liinsup{u{x) - v{y)) <F{p,Yo) - F{p,Xo) + K{C/{ay/^) 

where C is some constant and {p,Xo,Yo) is a limit point of {p,X,Y) as e | 0. If 
p ^ 0, we are done since 

F{P, Yo) - F{p, Xo) = F{p, Yo) - F{p, Xo) < 

because Xq < Yq. li p = 0, we use the information 4a|a; — y\'^{x — y) ^ p = 
and (9.7) (recall a is fixed) to conclude that Xq = Yq = 0, and then F{p,Yo) = 
F_{p, Xq) = 0, and we are still done. 

Proof that solutions lie in UC(R^). This is exactly as in §5.D. 

Proof of existence. A supersolution and subsolution are produced exactly as in 
§5.D. Perron's method still applies here, since F is lower semicontinuous and, we 
may use Lemma 4.4 as is with F replaced by F_ and the variant of Lemma 4.4 given 
in Remark 4.5 with G- = F and G+ = T. 

Notes on §9. As shown in §1, a number of equations arising from geometrical 
considerations present singularities at p = 0. The fact that this can easily be 
circumvented was shown independently by L. C. Evans and J. Spruck [72] and Y. 
Chen, Y. Giga, and S. Goto [45]. We also mention the work of H. M. Soner [152] 
on the equation (9.1) and the papers by G. Barles [16] and S. Osher and J. Sethian 
[140] , which showed how various geometrical questions about "moving fronts" could 
be reduced to equations that can be handled by viscosity theory. A general class 
of singular equations is treated in Y. Giga, S. Goto, H. Ishii, and M. H. Sato [80] 
and M. H. Sato [146]; these works establish existence and uniqueness as well as 
convexity properties of solutions. 

10. Applications and perspectives 

In this section we list some applications of the theory of viscosity solutions and 
indicate some of the promising directions for development of the theory in the next 
few years. We give some important references but they are not exhaustive. 

To begin this rather long list of applications, we recall that perhaps the main 
motivation for developing the theory was its relevance for the theories of Optimal 
Control and Diff'erential Games. Indeed, as is well known, in the theory of optimal 
control of ordinary differential equations or stochastic differential equations (with 
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complete observations) or in the theory of zero sum, two player deterministic or 
stochastic differential games, the Dynamic Programming Principle (DPP for short) 
states that the associated value functions should be characterized as the solutions 
of associated partial differential equations. These equations arc called Bellman or 
or Hamilton-Jacobi-Bellman (HJB for short) equations in control theory or Isaacs 
equations in differential games. The DPP was, however, heuristic and proofs of 
it required more regularity of the value functions than they usually enjoy. The 
flexibility of the theory of viscosity solutions has completely filled this regularity 
gap: roughly speaking, value functions are viscosity solutions and are uniquely 
determined by this fact (via the uniqueness of viscosity solutions). See, e.g., [70, 
116, 118, 119, 153]. This basic theoretical fact allows a spectacular simplification 
of the theory of deterministic differential games [70, 153] and also provided the 
possibility of creating sound mathematical foundations for stochastic differential 
games [77]. See also [88, 136] in addition to [29, 30]. 

The generalization of the definition of viscosity solutions to systems in diag- 
onal form is rather straightforward and has applications to optimal control and 
differential games. For these topics see, e.g., [43, 63, 92, 95, 114, 147]. In the 
case of systems, a combination of viscosity solutions and weak solutions based on 
distribution theory may define a natural notion of weak solutions [80] . 

More generally, as usual, a better understanding of existence-uniqueness issues 
for classes of equations leads to a better understanding of more specific issues. 
Typical examples here are perturbation questions, asymptotic problems, and a more 
detailed solution of some specific applications to Engineering or Finance problems 
[35, 161, 162]. Also, the part of the theory concerned with boundary conditions 
has led to a rather complete theory for problems with state-constraints (at least for 
deterministic problems) that are enforced by cost or boundary mechanisms [149, 
150, 44, 90, 25, 26, 33, 34]. Let us also mention that an interesting link between 
viscosity sohitions and the other main argument of Control Theory, the Pontryagin 
principle, has been shown [31, 9]. Of course, last but not least, these results have 
led to numerical approaches to Control or Differential Games problems via the 
resolution of the HJB (or Isaacs) equation. 

Indeed, the viscosity solutions theory is intimately connected with numerical 
analysis and scientific computing. First of all, it provides efficient tools to per- 
form convergence analyses (e.g., [52, 28, 133, 155]). It also indicates how to build 
discretization methods or schemes for other general boundary conditions and in 
particular for classical boundary conditions when working with rather degenerate 
equations [73, 144]. 

Another consequence of efficient existence, uniqueness, approximation, and con- 
vergence results is the possibility of establishing or discovering various qualitative 
properties of solutions (formulae, representations, singularities, geometrical prop- 
erties, characterizations and properties of semigroups, . . . — [12, 40, 41, 54, 56, 70, 
88, 124, 131]). Of course, one of the most important qualitative properties is the 
regularity of solutions. Viscosity solutions, because of their flexibility and their 
pointwise definition, have led to regularity results that are spectacular either in 
their generality (regularizing effects, Lipschitz regularity, or semicontinuity — [122, 
116, 57, 20, 18, 17]) or by their originality (C^'", C^'" or pointwise LP estimates— 
[39, 96, 158, 159, 160, 103]). 

The uniqueness and convergence parts of the theory have made possible partial 
differential equations approaches to various asymptotic problems like large devia- 
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tions [68, 76, 11, 141], geometrical optics [71, 23, 22, 13, 69], or homogenization 
problems [66, 123, 128] by arguments that are both powerful and simple. 

More specific applications concern the interplay between the behavior of so- 
lutions at infinity and structure conditions [83, 54, 56, 134] or the treatment of 
integrodifi^erential operators [112, 148]. Other applications concern some particular 
classes of equations arising in Engineering like some models of the propagation of 
fronts in Combustion Theory [16, 140], or the so-called shapes from shading models 
in Vision Theory [144]. 

Finally, a large part of the theory (but not yet all of it) has been "raised" 
to infinite-dimensional equations both for first-order and second-order equations 
[55, 125-127]. In addition to "standard" extensions to infinite-dimensional spaces, 
specific applications like the optimal control of partial differential equations or 
even stochastic partial differential equations — this last topic being motivated in 
particular by the optimal control of the so-called Zakai's operation, a well-known 
formulation of optimal control problems of finite-dimensional diffusion processing 
with partial observations — require some new developments of the viscosity solutions 
theory in order to accommodate unbounded terms in the equations [55, Parts IV 
and V; 126, 8, 7, 151, 42, 157]. 

The above rather vague and general comments on applications already contain 
many hints concerning promising directions the theory of viscosity solutions may 
take in the near future. In particular, the infinite-dimensional part of the theory 
will most probably explode in view of the unbounded possible avenues of inves- 
tigation. Much more progress is also to be expected for degenerate second-order 
equations and boundary conditions. In particular, progress is to be made on ex- 
istence questions for Dirichlet (and state constraints) boundary conditions — this 
might have a considerable impact on various applications like models in Finance. 
Similarly, we expect progress on uniqueness and regularity questions for uniformly 
elliptic second-order equations. 

It is also reasonable to hope that the insight gained by viscosity solutions will 
help to devise efficient high-order schemes for numerical approximations and prove 
their convergence. 

More specific developments should (and will) concern questions of behavior at 
infinity, geometrical optics problems, the use and the theory of discontinuous solu- 
tions, and the investigation of "second-order" integrodifferential operators .... 

Of course, the reader should not restrict his imagination to the borders we drew 
above. 

Appendix. The proof of Theorem 3.2 

In this section we sketch the proof of Theorem 3.2, which we reproduce here for 

convenience. 

Theorem 3.2. Let Oi be a locally compact subset of R^' for i = 1,. . . ,k, 

= 0ix---x0k, 

Ui e \JSC{Oi), and (p be twice continuously dijferentiable in a neighborhood of O. 
Set 

w{x) = ui{xi) -\ \-Uk{xk) for X = {xi, ■■■ ,Xk) e O, 
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and suppose x — (f i, . . . , Xk) E O is a local maximum of w — relative to O. Then 
for each £ > there exists Xi e S{Ni) such that 



where A = D^^{x) £ S{N), N = Ni + --- + Nk. 

The sketch we give below will be somewhat abbreviated, but the main points 
are all displayed and the reader should be able to provide any omitted details from 
the indications given. 

As defined and used in the main text, USC(O) consisted of the upper semicon- 
tinuous functions mapping O into M; however, at this point it is convenient to allow 
the value " — oo," so we hereafter set 

1180(0) = { upper semicontinuous functions u: O ^ RLi {—oo}} . 

The convenience is well illustrated by the following reduction. Below, conditions 
written for the index i arc understood to hold for « = 1, . . . , fc. 

First Reduction. We may as well assume that Oi — M^'. Indeed, if not, we 
first restrict Ui to a compact neighborhood Ki of Xi in Oi and then extend the 
restriction to M.^^ by Ui{xi) = — oo if ^ Ki (abusing notation and still writing 
Ui). The compactness of guarantees that Ui G USC(R^'); one then checks that 

2 + 2 + 

Jq Ui{xi) = J ' Ui{xi) (provided — oo < ^^(^j)). It is clear that x is still a local 
maximum oi w — Lp relative to M.^ = K^^^ This device of putting functions 

equal to — oo on unimportant sets will be used several times below to localize 
considerations to any compact neighborhood of a point under consideration. 

Second Reduction. We may as well assume that x = 0, Ui{0) = 0, (p{x) = 
^{Ax,x) for some A G S{N) is a pure quadratic, and is a global maximum of 
w — tp. Indeed, a translation puts x at the origin and then by replacing ip{x), Ui{xi) 
by (p{x) — (v(0) + {Dip{0),x)) and Ui{xi) — (ui(0) + (l?j,,<^(0), Xi)), we reduce to the 
situation x = D(p{0) = 0, and ip{0) = Ui{0) = 0. Then, since (p{x) = ^{Ax,x) + 
o{\x\'^) where A = D'^(p{0) and w{x) — (p{x) < w{0) — (p{0) = for small x, ifrj > 
we will have w{x) — ^{{A + r]I)x, x) < for small x ^ 0. Globality of the (strict) 
maximum at may be achieved by localizing via the first reduction. If the result 
holds in this case, we may use the bounds asserted in Theorem 3.2 to pass to the 
limit as ?7 J, to obtain the full result. 
Thus we have to prove the following 



Theorem 3.2'. Let m e USC(M^*). Ui{0) = 0, for i = 1,2, ... ,k, N = Ni + ■■■ + 

Nk, A e S{N), and 

(A.2) w{x) = ui{xi) H h Uk{xk) < ^{Ax,x) for x = {xi, ■ ■■ ,Xk) G M.'^. 



{DxiV{x),Xi) £ Jq Ui{xi) fori = l,... 



k 



and the block diagonal matrix with entries Xi satisfies 



(A.1) 
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Then for each e > 0, there exists Xi G S{Ni) such that 

2 + 

(0,Xj)€Jo Ui{xi) fori = l,...,k, 
and the block diagonal matrix with entries Xi satisfies 

'Xi ... 

(A.3) -ll + \\A\\]l< [ : •.. : I <A + eA^. 



£ 



... Xk. 



Here are the main steps in the proof: 
Step 1. Introduction of the sup convolutions. If e > 0, the Cauchy-Schwarz 
inequality yields 

(A.4) l^Ax,x)<{{A + EA^%Ci+(^^ + \\A\\^\x-e for,a;,eeK^. 
Putting 

(A.5) A = J + ||A|| 

and using (A. 2) and (A.5) we find 

(A.6) (^«i(a;i)-^ki-a|')+---+(M^fc)-^|a;fc-ai') <i((^ + £^')C,0 

or w(x) — (A/2)|a; — < (l/2)((^+£^^)^, ^). If necessary, we modify the definition 
of the Ui off a neighborhood of to assure that they are bounded above. Put 

(A.7) w{C) = sup (w{x) -^\x- , Uii^i) = sup (ui{xi) - ^\xi - 

xeRW \ / XiGR^i \ ^ 

SO that w{£_) = ui{£_i) + - ■ ■ + Uk{£.k)- Observe that, choosing Xi = 0, Ui{0) > Ui(0) = 
while ui(0) H h Ufc(O) < 0, so Wi(0) = 0. 

Step 2. Two results on semiconvex functions: Aleksandrov's Theorem 

and Jensen's Lemma. The functions iIi and Ui arc semiconvex; more precisely, 
since the supremum of convex functions is convex, w(^) + (A/2)|^p and Ui{£,) + 
(A/2)|^ip are convex. We then call A a semiconvexity constant for w and Ui. This 
is the reason we will employ two nontrivial facts about semiconvex functions. The 

first assertion is a classical result of Aleksandrov: 

Theorem A. 2. Let (p: R be semiconvex. Then cp is twice differentiable 

almost everywhere on M.^ . 

Proof. Since we assume that the members of our audience who are not functional 
analytically oriented have already seen a proof of this result, wc will give one for 
those who are so oriented. In the course of proof we assume certain facts about 
convex and Lipschitz functions that are known to most readers. References will be 
given in the notes. 
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Wc may as well assume that Lp is convex; a convex function on R is locally 
Lipschitz and hence once differentiable almost everywhere. We let 

Fi = {a; G : D^p{x) exists}; 

the letter F being used to indicate that F\ has full measure (i.e., its complement 
is a null set). The subdifferential d^p of y is the set- valued function given by 

y G d^{x) if ip{x) > (p{x) + {y,x — x) for x G . 

It is clear that d(p{x) = {D(p{x)} whenever ip is differentiable at x. In particular, 
d(p is single valued almost everywhere. By convex analysis, J = (I + dip)~^ is a 
single- valued mapping of into itself and is nonexpansive (has Lipschitz constant 
1). (Given x G R^, the solution 2; of 2; + d(f{z) 5 x is the value of J a.t x; z is 
the minimum of g — > (p{q) + ^\x — q\^.) Since (p is defined on all of M.^ , d(p{z) is 
nonempty for every z. We claim that for almost all x G M.^ Dip is differentiable at 
X, i.e., there exists A G S{N) such that 

Dip{y) = Dip{x) + A{y - x) + o{\y - x\) for y e 

for almost all x G M^. 
To establish this, let 

F2 = {J{x) : J is differentiable at x and DJ{x) is nonsingular}. 

Since J is Lipschitz continuous and onto (because d(p{z) is nonempty for every 
z), F2 has full measure. Wc use here the facts that {x : DJ{x) exists } is a set of 
full measure, J maps null sets to null sets and J{{x : DJ{x) exists and is singular}) 
is null. It is to be shown that D(p is differentiable on F3 = FiCi F2, which is a set 
of full measure. By the definition of J, 

Dip{J{x)) =x- J{x) for J{x) e F3C Fi. 

Thus assume J{x) + 5y £ Fi where 6y is small. By assumption, Dip{J{x) + 5y) 
exists. Since J is Lipschitz and DJ{x) is nonsingular, if 5y is sufficiently small, there 

is a 5x solving J{x + 5x) = J{x) + Sy, and wc may choose 5x to satisfy \Sx\ < K\5y\ 
for some constant K. Since J is a contraction, we also have \Sy\ < \dx\ and thus 
\dx\ and |i5!y| are comparable. Then 

D(p{J{x) + 6y) = D(p{J{x + Sx)) = x + Sx - J{x + 6x) 
= D(p{Jx) + {6x - DJ{x)Sx) + o{6x) 
= Dip{Jx) + {I - DJ{x))5x + o{5y). 

It remains to see that 5x — {DJ{x))~^5y + o{5y); however, this follows from the 
fact that \6x\ and \5y\ are comparable and the relation J{x) + DJ{x)Sx + o{Sx) = 
J{x) + Sy. We conclude that D{Dip){J{x)) exists and is DJ{x)~^ — I. 
It remains to establish that 



^{J{x) + Sy) = ^{J{x)) + Dip{J{x))5y + \{{DJ{x)-^ - I)5y, Sy) + oi\Syf) 
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for J{x) G F3. However, if we let tl'{5y) = (p{J{x) + Sy) and i'{Sy) = (p{J{x)) + 
D(p{Jx)Sy + ^{{DJ{x)~^ — I)Sy, 5y), we have V'(O) = ip{0) and for almost all small 
Sy 

Di,{5y) = Dip{J{x) + 5y) = D<fi{J{x)) + {DJ{x))-^5y - 6y + o{5y) 
= Di){5y) + o{5y); 

thus "^{Sy) =^ ip{Sy) — ijj{Sy) is locally Lipschitz continuous and satisfies ^'(0) = 
and D'^{5y) = o(5y) for almost all small 5y. It is clear that then '^{5y) = o{\5y\^), 
whence the result. 

The next result we will need concerning semiconvcx functions, which we call 
Jensen's lemma, follows. In the statement, B{x,r) is the closed ball of radius r 
centered at x and is the ball centered at the origin. 

Lemma A. 3. Let (p: R'^ M. be semiconvex and x he a strict local maximum 

point of ip. For p S M^, set ipp{x) = (p{x) + {p, x). Then for r, (5 > 0, 

K = \x € B{x,r): there exists p € Bg for which (pp has a local maximum at x} 

has positive measure. 

Proof. We assume that r is so small that (f has a; as a unique maximum point in 
B{x,r) and assume for the moment that if is C^. It follows from this that if S is 
sufficiently small and p € Bg, then every maximum of (pp with respect to B{x,r) 
lies in the interior of B(x,r). Since Dip + p = holds at maximum points of (/5p, 
Dip{K) D Bs- Let A > and ip{x)-^{X/2)\x\'^ be convex; we then have -A7 < D'^ip; 
moreover, on K, D^ip < and then 

-A/ < D'^ipix) < for X e K. 

In particular, | dct D'^(p{x)\ < for x G K. Thus 

meas{Bs) < meas{D(p{K)) < / | det D'^ ip{x)\dx < meas(Js:)|A|^ 

Jk 

(see the notes) and we have a lower bound on the measure of K depending only on 
A. 

In the general case, in which (p need not be smooth, we approximate it via 
mollification with smooth functions (p^ that have the same semiconvexity constant 
A and that converge uniformly to ip on B(x,r). The corresponding sets obey 
the above estimates for small e and 

00 00 

^ ^ n u ^1/™ 

n=l m=n 

is evident. The resiilt now follows. 

Step 3. A consequence of Step 2 and magic properties of sup convolution. 

Lemma A. 4 below applies to w of Step 1 and we shall see it provides us with matrices 
Xi G S{Ni) such that (0,Xj) G J^'"^{ti(0) and (A.3) holds. 
We use the notation 

(A.8) J\f{z) = n 

2 2 -|- 2 

from which one defines J analogously to J ' , J ' . Note that {p,X) G J'^f{x) 
amounts to 

fiy) = fix) + {p,y-x) + ^{X{y - x),y - x) + o{\x - yf) as y ^ x; 
i.e., / is twice differentiate at x and p = Df{x), X = D'^f{x). 
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Lemma A.4. /// e C(M^), B G S{N), /(0+(A/2)|^|2 convex anrf maxuiv (/(^)- 
1(5^,^)) = /(O), i/ien i/iere is an X e <S(iV) such that (0,X) e 7^(0) and 

-XI <X<B. 

Proof. Clearly f{£_) — ^{B£^,£^) — |^|^ has a strict maximum at ^ = 0. By Alcksan- 
drov's Theorem and the convexity assumption, / is twice difFerentiablc a.c. and 
then by Jensen's Lemma, for every S > there exists qs G with \qs\ < 6 such 
that /(^) + {qs,0 ~ ^{^^^0 ~ ICI^ has a maximum at a point with j^^l < S and 
is twice differentiable at ^g- Prom calculus, \qs\, \^s\ < S, and the convexity of 
/(0 + (A/2)|CP, 

Dfi^s) = 0{S), -\I < £>Vfe) <B + 0{5^). 

Observing that {Df{^s),D'^f{^s)) € J^f{xs), we conclude the proof upon passing 
to a subsequential limit (selected so that -0^/(^6) is convergent) as (5 | 0. 

Applying Lemma A. 4 to w from (A. 7) with B = A + eA^ and noting that the 
values of J^w) are formed as one expects given the representation w = u\-\ l-Wfe, 

2 

we obtain (0, X^) £ J ?ii(0), i — 1, . . . , /c such that (A. 3) holds. The final step in 
the proof of Theorem 3.2' is provided by applying the next lemma to each Wj to 
conclude that (0,Xi) e J '^Mi(O). 

Here are magical properties of the sup convolution. 

Lemma A. 5. Let A > 0, u € USC(IR^) hounded above, and 
v{C)= sup {v{x)-{\/2)\x-£,\^). 

//T^g e M^, y G S{M), and {q,Y) G J^'+v{r)), then 

iq,Y)GJ^'+v{rj + q/X) and vi?]) + il/2X)\q\'' = v{rj + q/X). 

2 . + 2 . + 

In consequence, if {0,Y) G J ' v{0), then (0, y) G J ' f(0). 

Proof. We assume {q,Y) G J^''^v{ri), and let y G be a point such that v{ri) = 
v{y) - (A/2)|y - rjf. Then for any x, ^ G 

v{x)-^\i-x\'<v{0 



(A.9) ^ ^ A 



= viy)-'-^\y-V? + {q,^-v) 

+ l{Y{i-v)i-v)+oi\^-V\') 
= v{y) - ^\y - + {q,^-v) + 0{\^ - vf). 
Putting ^ = x — y + r]m the appropriate one of these relations yields 
(A.IO) v{x) < v{y) + {q,x-y) + ^{Y{x -y),x-y)+ o{\x - y\^), 
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while substituting x = y and ^ = r] + a{X{r] — y) + q) yields 
(A.ll) 0<a\X{T]^y) + q\^ + 0{a^). 

The former simply says that {q, Y) e J'^'^v{y) while the latter with small a < 
implies that A(?7 — y) + g = or y = 77 + g/A as claimed. The relation v{r]) + 
(l/2A)|(7p = v{rj-\-q/X) follows at once. Assume now that F„) e J^'^w(^n) and 
{.U-,v{i),qn,yn) (0,0, 0,y); by the foregoing, (g„,F„) £ J^'+w(C„ + Qn/A) and 
v{£,n + </n/A) = t'(Cn) + (1/2A) ^ . By the definitions, v{0) < v{0) and from this, 
the upper semicontinuity and the foregoing, we have 

v{0) > limsupz;(^„ + ^) = limsup(6(e„) + ^\qn\^) = ^(0) > v{0), 

which provides the final piece of information we needed to conclude that (0, F) G 

Remark A. 6. We close with a final and more elegant reformulation of what was 

proved above: if the Ui satisfy the conditions of Theorem 3.2, w{x) = Ui{xi) + 

2 + 

• • • + Uk{xk) and {{pi, ■ ■ ■ ,pfe), A) e J^' ^^(a^)) then for each £ > 0, there exists 

2 -|- 

Xi G S{Ni) such that e Ui{xi) and (A.l) holds. 

Notes on the appendix. Above, we provided (for the first time) a self-contained 

proof of Theorem 3.2. Except for the two auxiliary results on semiconvex functions, 
the main tool is the so-called sup convolution. This approximation procedure (more 
often in the guise of inf convolution) is well known in functional analysis and, in 
particular, in convex analysis and the theory of maximal monotone operators (see, 
for example, the text [38] of H. Brezis). It was noticed in [110] that it may provide 
an efficient regularization procedure for (even degenerate) elliptic equations; some 
of its properties are given there. See also [104]. Its "magical properties" can be 
seen as related to the Lax formula for the solution of 

du) X 

— - -|Vw|^ =0 for a; e K^, i > 0, w\t=o = v on , 
at 2 

which is 

w{x,t) = sup|t;(y) - ^l^^-ypj • 

Indeed, the coincidence of this solution formula and solutions produced by the 
method of charteristics leads to the properties used. Of course, this is a heuristic 
connection, since characteristic methods require too much regularity to be rigorous 
here. 

The inf convolution can also be seen as a nonlinear analogue of the standard 
mollification when replacing the "linear structure of and its duality" by the 
"nonlinear structure of or C." One can also interpret this analogy in terms of 
the so-called exotic algebra (K, max, -|-). 

Theorem A. 2 is a classical result of A. D. Aleksandrov [1]. The proof given here 
is in fact a slightly stronger form of this result {D(f is difi'erentiable a.e.) and our 
proof follows F. Mignot [135] (although we have made his proof more complicated 
for pedagogical reasons). In the proof we used the fact that Lipschitz functions are 
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difFerentiable a.e., which is called Rademacher's Theorem. A proof may be found in 
L. C. Evans and R. Gariepy [67] and F. Mignot [135]. We also used that Lipschitz 
functions map null sets to null sets and the solvability of J{x + 5y) = J{x) + 5x 
when DJ{x) exists and is nonsingular. These are proved in, respectively, Lemma 
7.25 and the proof of Theorem 7.24 in W. Rudin [145]; see also [135]. We also 
used the fact that J{{x : DJ{x) exists and is singul^) is null. This follows from 
the general formula 

/ #{AnJ-^{y))dy= [ \ det DJ{x)\dx 

where # is counting measure. This formula holds for Lipschitz continuous func- 
tions and measurable sets A C and is a special case of the Area Formula for 
Lipschitzian maps between Euclidean spaces. The Area Formula may be found in 
L. C. Evans and R. Gariepy [67] and H. Federer [74]. 

Jensen's Lemma (Lemma A. 3) is a variation on an aspect of the theme known 
as "Aleksandrov's maximum principle" (see [2, 3, 6, 37, 120, 143]). 
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